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SUMMARY 


I 


/ 

S 

'  The  problem  of  the  active  control  of  linear  elastic  structures 
using  piezoceramic  transducers  as  sensors  and  actuators  has  been 
investigated  by  a  combined  theoretical  and  experimental  approach.  The 
optimal  rate  feedback  gain  distribution  of  an  active  structure  with 
multiple  collocated  sensors  and  actuators  has  been  obtained  by  using  a 
limited  state  feedback  approach  which  resulted  in  an  increase  in 
system  damping.  To  model  the  active  structure  for  the  optimal  control 
problem,  a  finite  element  model  has  been  developed.  An  active  element 
consisting  of  a  simple  beam  element  with  a  bonded  unimorphic 
piezoceramic  sensors  and  actuators  has  been  obtained.  The  model 
incorporates  the  electromechanical  coupling  of  the  transducers, 
bonding  effects  and  a  mathematical  model  for  the  feedback  signal 
conditioning  circuitry.  The  resulting  discrete  degrees  of  freedom 


model  is  in  the  form  of  a  set  of  coupled  ordinary  differential 
equations  which  describe  the  dynamic  behavior  of  the  active  structure. 
To  obtain  the  unknown  dynamic  coupling  coefficients  that  represent  the 
effects  of  bonding  and  other  parameters  of  the  model  accurately, 
parameter  identification  methods  have  been  used.  ^Modal  control  has 
also  been  experimentally  demonstrated  by  conditioning  the  output  of 
each  individual  sensor  with  an  adjustable  bandpass  filter,  phase 
shifting  and  gain  circuits.  The  identification  of  bi orthogonal  modes 


of  the  resulting  non  self  adjoint  system  when  non  collocated  sensors 
and  actuators  has  also  been  accomplished.  The  identified  discrete 


••.'a'* 


degrees  of  freedom  model  and  a  quadratic  performance  Index  have  been 
used  in  obtaining  optimal  control  laws.  In  this  phase,  the  problem 
has  been  treated  as  a  regulator  with  limited  state  feedback.  As  a 
next  step,  the  optimal  control  problem  has  been  solved  by  considering 
the  active  structure  as  a  distributed  parameter  system.  An  optimal 
control  law  has  been  obtained  by  maximizing  the  decrease  of  the  time 
derivative  of  a  Lyapunov  functional  of  a  cantilever  beam  with  a 
collocated  sensor  and  actuator  occupying  a  subdomain  of  the  structure 
with  rate  feedback  control.  The  developed  control  law  has  been 
validated  by  using  an  explicit  finite  difference  method.  The 
governing  partial  differential  equations  have  been  solved  for  the 
system  subject  to  excitations  and  control. 
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CHAPTER  I 

INTRODUCTION  AND  BACKGROUND 

In  the  past  few  years,  there  has  been  considerable  research 

activity  in  the  field  of  active  and  passive  control  of  vibrations  of 

flexible  structures.  One  of  the  methods  of  active  control  of 

vibrations  involves  the  placement  of  piezoelectric  devices  on  a 

structure  to  sense  and  control  dynamic  strains  induced  by  structural 

vibrations.  The  deformation  of  the  sensing  transducer  results  in  an 

electric  current  which  is  conditioned  by  operations  such  as 

amplification  and  shifting  of  the  phase  of  the  signal.  The  conditioned 

signal  is  then  applied  to  another  piezoelectric,  electrostrictive  or 

magnetostrictive  device  placed  at  a  selected  location  of  the 

structure.  This  transducer  acts  as  an  actuator  motor  or  a  driver  and 

transmits  mechanical  energy  to  the  structure.  Depending  on  the 

applied  voltage,  electromechanical  coupling  of  the  forcing  transducer 

to  the  structure  and  its  location,  a  certain  amount  of  structural 

control  of  flexible  structures  can  be  achieved.  The  term  "electronic 

damping  "  has  been  used  to  describe  this  method  of  structural  control 

1-5 

in  some  of  the  more  recent  literature. 

This  type  of  active  control  offers  unique  features  that  are  not 
usually  employed  for  control  of  structural  vibrations.  For  example, 
since  the  dynamics  of  direct  contact  type  sensors  and  actuators  permit 
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a  wide  frequency  range  of  control,  a  measure  of  tunability  is  provided 

for  the  control  of  structural  systems  that  age  or  grow.  This  method 

also  has  the  potential  benefit  of  adding  less  mass  to  the  system  than 

comparable  passive  and  active  servo  controller  systems. 

Passive  systems  lack  narrow  band  tunability  and  active  servo 

controllers  require  an  inertial  ground  to  react  against.  It  is 

expected  though  that  each  particular  type  of  vibration  control 

device  will  be  beneficial  in  different  areas  of  operation  in  terms  of 

the  amplitude,  control  force  requirements  and  frequency  bandwidth. 

The  potential  of  using  piezoceramic  transducers  for  vibration 

control  has  been  convincingly  demonstrated  in  a  variety  of 

2  5 

experiments  such  as  jitter  control  of  optical  systems.  ’ 

Nevertheless,  there  are  many  issues  which  still  need  to  be  resolved 
before  appropriate  control  design  procedures  are  developed. 
Specifically,  given  a  piezoelectric  material,  electro-mechanical 
coupling,  the  structure  and  the  control  domain,  a  designer  should  be 
able  to  obtain  an  optimum  procedure  for  controlling  the  vibration  of 
the  given  system.  To  meet  this  design  goal,  one  must  be  able  to  obtain 
a  measure  of  the  system  changes  in  damping  and  stiffness  due  to  the 
active  feedback  control  similar  to  the  way  the  mass  and  stiffness 
properties  are  obtained  by  using  other  analytical  techniques. 

To  resolve  some  of  these  issues,  the  research  in  this  thesis  is 
directed  towards  the  following  major  objectives; 

(1)  The  first  objective  is  to  develop  models  for  the  active  structure 
consisting  of  a  beam  structure,  piezoceramic  sensors,  signal 


conditioning  systems  with  appropriate  control  laws  and 
piezoceramic  actuators.  The  term  “active  beam"  or  “active 
structure"  has  been  used  to  refer  to  such  a  structure. 

(2)  The  second  objective  is  to  identify  the  dynamic  characteristics 
of  the  active  structure  from  measurements  of  force  and 
acceleration  of  the  active  structure.  Cases  of  collocated  and 
non  collocated  sensors  and  actuators  will  be  considered. 

(3)  The  third  objective  is  to  identify  the  system  parameters  from 
measurements  of  the  input  and  output  of  the  active  structure.  In 
particular,  electromechanical  dynamic  coupling  coefficients  and 
structural  damping,  that  cannot  be  determined  otherwise,  are 
identified. 

(4)  The  final  objective  is  to  develop  optimal  control  algorithms  to 
selected  cases  involving  an  active  beam  with  piezoceramic  sensors 
and  actuators. 

In  the  next  section,  a  brief  review  of  the  field  of 
piezoelectricity  is  described. 

1.1.  Background 

Piezoelectricity  was  first  discovered  by  the  Curie  brothers.® 
This  work  led  them  to  look  at  the  electrification  of  crystals  upon  the 
application  of  pressure.  The  converse  piezoelectric  effect,  which  is 
the  change  in  crystal  dimensions  upon  the  application  of  an  electric 
field  was  first  theoretically  predicted  by  Lippman  and  in  the  same 
year  verified  by  the  Curies.®  Later  investigators  established  the 
relationship  between  piezoelectricity  and  the  crystalline  structure  of 


various  materials. 


The  first  practical  application  of  the  piezoelectric  effect  was 
for  the  development  of  acoustic  waves  in  fluids.  Quartz  crystalline 
material  was  initially  used  until  the  development  of  the  ferroelectric 
class  of  materials.®  The  term  ferroelectric  is'  used  to  describe  a 
material  which  does  not  have  a  center  of  symmetry  or  electric  dipole 
moment  above  a  certain  temperature  called  the  Curie  temperature.^ 

The  single  most  important  factor  of  the  ferroelectric's  over 
other  piezoelectric  materials  is  that  they  have  a  higher 
electromechanical  coupling.  Electromechanical  coupling  is  defined  as 
the  amount  of  deformation  achieved  per  field  strength  applied  to  the 
piezoelectric  material. 

Ferroelectric  materials  contain  both  organic  and  inorganic 
compounds.  The  piezoceramic  ferroelectric  materials  are  used  in  this 
investigation  for  structural  control  applications. 

The  constitutive  relationships  for  piezoceramic  material  have 
been  discussed  by  a  number  of  investigators  over  the  years  and 
compiled  by  the  IEEE  and  the  Electronic  Industries  and  American 

Q  Q 

Standards  Association  into  standards.  ’  One  of  the  most  serious 
disadvantages  of  piezoceramic  ferroelectric  type  materials  with 
reference  to  this  application  is  the  sensitivity  of  their 
piezoelectric  properties  to  changes  in  temperature.  D.  Berl incourt^® 
has  discussed  the  relationship  of  the  performance  of  selected 
piezoceramics  with  temperature  and  other  environmental  factors. 
Doping  of  piezoceramics  with  lead  and  other  elements  helps  to 
stabilize  these  parameters  over  a  wider  temperature  range.  The  lead 
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zirconate  titanate  (PZT)  type  piezoceratmcs  have  proven  to  be  the  most 
promising  compounds  for  stable  piezoelectric  properties  and  mechanical 
behavior  over  a  temperature  range  that  Includes  -200  to  150  degrees  C. 

Organic  type  piezoelectric  compounds  such  as  polyvinyl  1  dene 

fluoride  film  (PVDF)^^  have  been  used  for  active  control  of 

structures.  This  type  of  organic  ferroelectric  material  has  the 

potential  of  providing  strains  of  an  order  of  magnitude  greater  than 

piezoceramics  but  at  the  cost  of  a  much  greater  applied  electric 

field.  Some  of  the  more  recent  applications  of  PVDF  has  been  in  the 
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area  of  adaptive  optics  ’  and  in  vibration  control  of 
14  15 

structures."^  ’ 

The  electrostrictive  or  magnetostrictive  effects  in  some 
materials  can  perform  the  same  function  as  the  converse  piezoelectric 
effect.  Electrostriction  is  related  to  a  material  change  of  shape  upon 
the  application  of  a  large  electric  field.  Usually  the  change  in  shape 
is  very  small  and  until  recently  was  considered  as  a  second  order 
effect.  W.P.  Mason^®  has  studied  the  electrostrictive  effect  in 
nonpolarized  ferroelectric  type  ceramic  of  materials.  According  to 
Mason  the  electrostrictive  effect  is  quite  large  when  compared  to  the 
corresponding  converse  piezoelectric  effect  of  the  same  polarized 
material.  One  of  the  advantages  of  electrostrictive  actuators  is  that 
they  do  not  exhibit  hysteresis  effect  common  to  piezoceramic  type 
materials. This  feature  may  make  electrostrictive  ceramics  to  be 
more  suitable  as  driving  devices  for  vibration  control. 
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Some  of  the  applications  of  piezoceramics  can  be  classified  as 

quasistatic  applications,  dynamic  applications  with  resonanting 

piezoceramics  and  non  resonanting  dynamic  applications.^  Most  of  the 

work  involving  quasistatic  applications  is  in  the  area  of  adaptive 

optics.  Dynamic  applications  with  resonanting  piezoceramics 

have  been  primarily  used  in  creating  acoustic  surface  and  bulk 

waves.  Non  resonating  dynamic  application  are  mostly  in  the  area 

of  transducers  and  only  recently,  non  resonating  dynamic  applications 

have  been  in  the  area  of  active  control  of  flexible  structures. 
1-5,15,32,33 

Electronic  Damping 

One  of  the  earliest  known  uses  of  piezoceramics  to  active  control 
of  structures  is  attributed  to  Olsen. ^  In  this  work  he  proposed  an 
electronic  vibration  reducer  consisting  of  a  piezoelectric  driver  and 
sensor  attached  to  a  structure  with  a  suitable  feedback  amplifier. 

32 

Another  early  example  of  this  type  of  control  is  by  Mckechnie 

who  has  used  piezoceramic  sensors  in  a  specially  built  accelerometer 

to  drive  a  forcing  piezoceramic  to  reduce  the  resonant  response  of  an 

accelerometer.  His  experimental  results  indicate  a  reduction  in  the 

predicted  peak  response  by  thirty  percent. 

The  first  application  to  the  active  control  of  a  structure  is  due 
2 

to  Forward  who  has  investigated  the  vibration  control  of  a  mirror 


subjected  to  acoustical  excitation  of  its  resonances.  One 
piezoelectric  transducer  has  been  used  to  sense  the  vibrations  and  its 
output  has  been  conditioned  by  a  negative  feedback  amplifier  to  drive 


another  spatially  separated  transducer.  The  vibration  response  of  a 
single  mode  of  the  optical  structure  was  reduced  significantly. 

3 

Forward  and  Lui  have  also  demonstrated  additional  applications 
of  "electronic  damping"  to  control  resonances  in  gimbal  torquer 
control  loops  where  they  reduced  the  peak  response  of  a  gimbal  torquer 
loop  at  various  modal  frequencies  corresponding  to  gimbal  mechanical 
resonances.  This  was  accomplished  by  using  small  piezoelectric  strain 
transducers  to  provide  signals  into  a  summing  junction  in  the  torquer 

4 

control  loops.  Forward  and  Swigert  have  applied  an  "electronic 

damping"  system  to  reduce  orthogonal  bending  modes  in  a  cylindrical 

satellite  antenna  mast.  The  application  of  piezoceramic  sensors  and 

actuators  to  reduce  the  amplitude  of  response  of  two  problem  modes  in 

a  large  composite  optical  bench  has  been  reported  in  reference.^  One 

piezoelectric  ceramic  strain  transducer  has  been  used  as  a  sensor,  and 

ten  others  were  wired  in  parallel  as  a  combined  driver  in  a  velocity 
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feedback  loop.  In  a  recent  work,  Crawley,  has  presented  a 
mechanics  coupling  type  of  model  for  a  piezoceramic  bonded  to  a  small 
cantilever  beam.  The  model  predicts  the  first  mode  dynamic  response 
of  a  beam  when  driven  by  the  piezoceramic  transducers. 

An  examination  of  the  reported  literature  in  the  field  of  active 
control  of  structures  with  piezoceramic  sensors  and  actuators  has 
revealed  the  following.  A  consistent  analytical  model  does  not  exist 
for  an  active  structure  with  piezoceramic  or  piezoelectric  sensors  and 
actuators  that  occupy  selected  subdomain  of  the  structure.  The 
problem  of  identification  of  the  dynamic  characteristics  of  an  active 


structure  from  measurements  of  time  histories  of  appropriate  input  and 

output  has  not  been  explored.  In  particular,  a  significant  problem  of 

practical  importance  is  the  identification  of  dynamic  characteristics 

of  an  active  structure  with  non  collocated  sensors  and  actuators. 

Procedures  for  the  identification  of  system  parameters  such  as 

electromechanical  coupling  does  not  exist.  The  problem  of  optimal 

control  of  structures  with  piezoceramic  sensors  and  actuators  has  not 

been  examined.  Only  the  problem  of  “instantly  optimal control  of  a 

structure  with  a  single  PVDF  actuator  and  a  conventional  accelerometer 

15 

sensor  has  been  explored. 

Some  of  these  unsolved  problems  discussed  in  the  previous 

paragraph,  have  been  examined  in  the  subsequent  chapters.  In  chapter 
II,  the  development  of  a  finite  IV  a  method  is  developed  for  model 
identification  for  a  structure  with  non  collocated  sensors  and 
actuators.  The  optimal  gain  distribution  of  multiple  location  control 
of  a  simple  structure  using  rate  feedback  is  explored  by  applying 
limited  state  output  feedback  control  in  chapter  V.  Finally,  the  time 

optimal  control  of  a  cantilever  beam  using  a  distributed  piezoceramic 

sensor  and  actuator  is  presented  in  chapter  VI. 


CHAPTER  II 


FINITE  ELEMENT  MODEL  OF  THE  ACTIVE  STRUCTURE 

In  this  chapter  a  finite  element  model  has  been  developed  for  an 
active  structure  shown  in  figure  2.1  consisting  of  a  linear  isotropic 
elastic  beam,  piezoceramic  sensor  occupying  a  subdomain  of  the  beam, 
with  the  signal  conditioning  and  piezoceramic  actuator  occupying  another 
subdomain  of  the  beam.  In  the  following  sections,  the  sensor,  actuator 
and  control  electronic  are  analyzed  and  the  finite  element  model  for  the 

active  structure  has  been  developed.  Previous  studies  have  considered 
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only  a  finite  element  model  or  a  beam  with  a  piezoelectric  actuator'^^’  . 

2.1.  Sensing  Transducer 

The  sensing  transducer  consisting  of  a  piezoceramic  thin  plates  is 
assumed  to  be  bonded  to  a  beam  structure.  The  basic  construction  of  the 
transducer  consists  of  a  thin  layer  of  G-1195  PZT  ceramic  coated  with  a 
thin  nickel  electrode  film.  In  this  form  the  piezoceramic  transducer  is 
called  a  unimorph.  When  two  unimorphs  are  bonded  together  and  share  the 
common  electrode  between  them  they  are  called  bimorphs.  The 
polarization  of  the  unimorph  or  bimorph  depend  on  the  intended 
application. 

A  idealized  sketch  of  a  section  of  a  beam  with  a  sensing  transducer 


and  signal 

conditioning 

circuit  is 

shown 

in 

Figure  2.2. 

The 

piezoceramic 

transducer 

is  polarized 

in 

the 

+z  direction 

and 

electrically  connected  via  two  short  leads  a  current  amplifier.  The 
following  assumptions  are  made: 


a)  Initially,  a  perfect  bond  is  assumed  between  the  elastic  body 
and  the  transducer.  In  reality  there  is  a  relative  displacement  due  to 
the  deformation  of  the  bonding  material.  This  has  been  initially 
neglected.  Later,  this  effect  is  included  approximately  in  the  form  of 
a  dynamic  coupling  coefficient. 

b)  The  thin  film  electrode  surfaces  add  no  mass  or  stiffness  to  the 
transducer  and  the  transducer  itself  does  not  add  any  mass  or  stiffness 
to  the  structure  to  which  it  is  bonded. 

c)  The  capacitance  of  the  transducer  leads  is  considered  to  be 
negl igible. 

Considering  a  one  dimensional  analysis  of  the  sensor,  the  following 
constitutive  relations  derived  in  Appendix  A  are  approximated  for 


Sjj{x,t)  =  SjjTjj(x,t)  +  dj^EjCx.t)  (2.1) 

D3(x,t)  =  d3jTjj(x,t)  +  e[3E3(x,t)  (2.2) 


where 

^11  "  strain 
^11  “  stress 

E3  =  electric  field  between  the  electrodes 


D3  =  charge  per  unit  area 
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*11 

^31 

T 

=•13 


elastic  compliance  at  constant  electric  field 

piezoelectric  constant 

dielectric  constant  at  constant  stress 


Notations  used  in  eqs.  (2.1)  and  (2.2)  are  based  on  the  standards 
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developed  for  piezoelectricity  *  \ 

The  sensor  thickness  is  an  order  of  magnitude  smaller  then  the 
elastic  body  and  the  strain  distribution  is  assumed  to  be  constant 
through  the  sensor  cross  section  and  equal  to  the  upper  fiber  strain  (@ 
z=h)  of  the  elastic  beam.  Therefore  the  strain  distribution  along  the 
sensor  is  approximated  as 


Sij(x.t)  =  -  h 


(2.3) 


3X‘ 


When  the  sensor  is  deformed,  both  a  charge  and  an  electric  field  is 
produced.  The  electric  field  between  the  two  electric  plates  results  in 
an  electric  potential  at  any  cross  section  equal  to 


Vs(Xi,t)  =  J  E(x.,t)  dz 


(2.4) 


or 


Vg(x.,  t)  =  6  E(x^,  t)  (2.5) 

where  6  is  the  thickness  of  the  sensor.  Since  there  is  no  input 

go 

resistance  to  node  1  in  Figure  2.2  the  following  conditions  arise  , 


'‘A 

% 


vs 


•v' 


Cv*' 


f'*i 


|v^{t)  -  v_(t)|  = 


(2.7) 


for  a  typical  commercial  operational  amplifier  c^is  of  the  order  of  10  ^ 
volt.  Therefore  the  potential  difference  between  the  sensor  electrodes 


Vg(t)  =  1.0  X  lo‘\olt. 


(2.8) 


For  the  transducers  used,  6  is  approximately  2.54x10'^  M  therefore  the 
field  between  the  plates  is  of  the  order  of 


E3(t)  =  4.0  V/M  . 


(2.9) 


Multiplying  the  field  by  d^j  and  in  eqs.  (2.3)  and  (2.4)  a 
strain  of  7.16  x  10  M/M  and  electric  displacement  of  6.6  x  lO"®  C/M*M 
are  obtained.  These  terms  are  negligible  compared  to  the  strain  and 
charge  produced  due  to  the  stress  in  the  piezoceramic.  Furthermore 
because  of  the  negligible  electric  field,  the  sensor  ideally  becomes  a 
pure  current  source.  With  this  particular  configuration  using  eq. 
(2.3)  and  solving  for  the  stress  in  eq.  (2.1)  yields 


j  (x  t)  “ - — 

11'*’  —  S  E  2 


(2.10) 


and  the  electric  displacement  in  eq.(2.2)  as 


(2.11) 


Hence  the  charge  developed  related  to  the  sensor  strain  at  every 
point  given  a  constant  width  transducer  is 


•  A 


DjCx.t)  =  -  hd2j  Cjj 


E  a^wfx, 


(2.12) 


where  c^j  =  l/s^^^is  now  the  elastic  modulus  of  the  piezoceramic  material 
measured  at  constant  field  conditions.  The  input  sensor  current  ig{t)  is 
dependent  on  the  rate  of  charge  developed 

digCt)  =  fi3(t)dxdy.  (2.13) 
For  a  constant  width  plate  the  current  is 


.(t)  =  J  63  dx 


(2.14) 


where  w^  is  the  width  and  L  is  the  plate  length. 

The  charge  D3  in  general,  is  a  function  of  x  and  time.  Eqs.  (2.12) 
and  (2.14)  are  combined  to  yield 


i  =  -w  hd  c^  f  (x,t)  j 

^s  ”d"°3l‘^ll  J  j^2 


Using  Kirchhoff's  current  law  for  node  1  in  Figure  2.2 


(2.15) 


is  =  V  ^  i. 


furthermore,  for  the  circuit  shown  in  Figure  2.2 

M^it)  =  -Rpip(t) 

and  since  for  an  ideal  op  amp  (34) 


(2.16) 

r;’*i 

.s' ' 

(2.17) 
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i  _  =  0 

the  output  voltage  \/o(t)  becomes 


(2.18) 


VgCt)  =  -  Rpi3(t).  (2.19) 

The  final  result  is  obtained  by  combining  eqs.  (2.17)  and  (2.21) 


Vo(t)  = 


r  ’  ^  (x, 

J  3x2 


t) 


dx 


(2.20) 


M-1 


where 


Equation  (2.22)  represents  the  electromechanical  equation  for  the 
sensor  perfectly  bonded  to  an  elastic  beam.  It  provides  the  voltage 
output  of  the  sensor  as  a  function  of  the  strain  distribution. 
Limitation  of  the  Distributed  Sensor  Model 

Consider  a  unimorph  piezoceramic  sensor  bonded  along  the  entire 
length  of  a  simply  supported  beam  as  shown  in  Figure  2.3  where  h  is  the 
distance  from  the  neutral  axis  to  the  upper  beam  surface.  The  motion  of 
the  beam  under  conditions  of  free  vibrations  is  approximated  by 


w  (x,t)  "  ^  sin  sin  n^t  . 


(2.21) 


n=l 


By  considering  only  the  first  term  of  the  series 


v.-\ 
-vO 
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From  eqs.  (2.20)  and  (2.22) 

n=l;  V^(t)  =  2Kg  f  cos  t  .  (2.23) 

For  the  second  term  of  the  series  gives 

n  =  2:  \l^(t)  =  0  .  (2.24) 

Under  pure  second  mode  excitation,  ideally  half  of  the  sensor 
would  be  under  compression  and  half  of  it  would  be  under  tension 

resulting  in  a  zero  output  voltage.  Therefore  the  limitations  of  the 
sensing  model  is  that  the  highest  frequency  for  control  is  limited  to 
that  below  the  modal  density  which  results  in  a  nodal  line  within  the 
domain  of  the  sensor. 

2.2.  Actuator  Transducer 
35 

In  a  previous  study  ,  a  finite  element  model  for  a  beam  with  an 
electrostrictive  transducer  has  been  derived  using  an  energy  approach. 
The  electrostrictive  transducer  behaves  similarly  to  the  piezoelectric 
transducer  with  a  driver  with  an  applied  voltage.  In  this  thesis  a 

model  is  derived  using  a  force  equilibrium  approach  (Appendix  B)  and 
agrees  with  the  previous  work. 

The  piezoceramic  actuator  configuration  is  basically  the  same  as  a 
sensor  except  that  in  this  case  a  relatively  high  voltage  of  the  order 
of  100  to  120  volt(RMS)  can  be  applied  to  the  piezoceramic  plates.  When 
a  high  voltage  is  applied,  the  converse  piezoelectric  effect  occurs 

which  results  in  a  deformation  of  the  piezoceramic  material.  A 
particular  unimorph  can  be  polarized  such  that  a  positive  electric  field 
results  in  an  elongation  in  the  directions  perpendicular  to  the 
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polarization  vector  with  a  simultaneous  contraction  in  directions 
parallel  to  the  polarization  vector.  A  negative  voltage  would  result  in 
the  opposite  responses. 

When  an  actuator  transducer  is  bonded  to  an  elastic  structure  it  is 
constrained  and  as  an  electric  field  is  applied  the  actuator  tends  to 
elongate  or  contract  the  resulting  body. 

A  composite  structure  consisting  of  an  elastic  beam  and  two 
piezoceramic  actuators  is  shown  in  figure  2.4.  The  two  unimorph 
actuators  bonded  to  the  top  and  bottom  of  a  thin  beam  and  wired  in  such 
a  manner  that  when  a  voltage  is  applied  to  the  transducers  the  bottom 
actuator  expands  and  the  top  actuator  contracts.  To  analyze  this 
sandwich  construction,  an  inverted  form  of  the  constitutive  equations 


2.1  and  2.2  are  (Appendix  A) 


Tjj(x,t)  =  CjjSjj(x,t)  -  H3jD3(x,t) 


E3(x,t)  =  -  H3jSjj(x,t)  +  D3(x,t) 


(2.25) 


(2.26) 


where 


=  Strain 
=  Stress 

=  Charge  per  unit  area 

=  Elastic  modulus  with  the  piezoceramic  circuit  open 
=  Dielectric  constant  at  constant  strain 


I 

J 


=  Electric  Field  . 

Since  the  drivers  are  fully  covered  with  electrodes  the  electric  field 
E^  produced  between  the  plates  for  a  given  potential  applied  can  be 
considered  only  dependent  on  time. 

Because  of  the  direction  of  polarization  E^  is  assumed  to  be  much 
greater  than  Ej  and  E2.  Since  E^  »  Ej  and  E^  »  E2,  assuming  E^  =  E2  = 


(2.27) 


Vp(t)  =  E3{t)6 


where  6  is  the  separation  between  the  electrode  plates. 


(2.28) 


The  potential  applied  across  the  plates  of  the  top  unimorph  is  then 


equal  to 


II  » II  *  w 

=  j  E3(t)dz  =  I  (-H31SJJ  +  l/E33D3)dz  .  (2.29) 


Furthermore  when  pure  bending  is  assumed. 


S,,(x,t)  =  -  z  ^ 


(2.30) 


and  the  piezoceramic  is  a  dielectric  material  that  is  polarized  in  the  z 
direction,  the  charge  within  the  body  is  governed  by  (Appendix  A) 


Using  conditions  given  by  eqs.  (2.30)  and  (2.31) 
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Given  two  unimorph  actuators  wired  as  previously  mentioned,  the 
moment  for  at  any  cross  section  of  the  beam  of  unit  thickness  is  given 


M(x,t)  =-  J  Ez^  dz  +  j  dz  +  I  ^u^dz 
-h  -{h+6)  h 


(2.37) 


ii 

•j. 


r 

*  •  ^  l«  •  % 


2..  e-,z 


M{x,t)  =  El  ^  -  2  I  (cJjZ^  \/^(t))dz  .  (2.38) 


Neglecting  terms  that  are  multiplied  by  6  and  6*6  the  final  result  is 


^  ..  ^ 


M(x.t)  =  El  ^  253^h  Vp(t) 

dX 

where  e^j  is  equal  to  e^^  times  the  width  of  the  actuator, 

If  only  one  surface  has  an  actuator  the  moment  equation  is 

2 

M(x.t)  =  EI^  KMt) 

3x 


where 


(2.39) 


(2.40) 


r"  1 


Kq  = 


From  eq.  (2.40)  certain  observations  are  made.  As  a  voltage  is 
applied  to  the  piezoceramic  unimorph  it  is  equivalent  to  adding  opposing 
couples  at  the  ends  of  the  section  whose  sign  is  dependent  on  the  sign 
of  the  applied  potential  to  the  plates^^’^®. 


'■  1. 


Consider  a  structure  as  shown  in  figure  2.5,  with  no  external 


forces.  For  an  applied  DC  voltage  to  the  actuators  eq.  (2.40)  becomes 


^d'^d 


El 


3^w 

3X^ 


(2.41) 


Integrating  eq.  (2.41)  twice  and  solving  for  the  unknown  constants  using 
the  boundary  conditions  at  x  =  0  and  x  =  +  L  gives 


w  = 


El 


(2.42) 


For  the  same  structure  without  unimorph  piezoceramic  layers  and  loaded 
as  shown  in  figure  2.6,  the  resulting  expression  for  displacement  is 


w  = 


(2.43) 


A  comparison  eqs.  (2.42)  and  (2.43)  further  supports  the  observation 
that  a  piezoceramic  unimorph  driver  produces  equivalent  point  or 
concentrated  external  moments  proportional  to  the  voltage  applied. 

2.3.  Analysis  of  Active  Electronic  Feedback  Damping  Using 
the  Finite  Element  Method 

A  finite  element  approximation  of  the  coupled  sensor,  actuator  and 
structure  is  developed  in  this  section.  In  previous  studies,  finite 
element  models  for  such  an  active  structural  dynamic  system  have  not 
been  developed.  To  accomplish  this  a  beam  with  only  a  piezoceramic 
unimorph  actuator  is  considered  first.  Later  an  element  consisting  of  a 
collocated  sensor  and  actuator  transducer  is  examined. 


The  beam  drawn  in  figure  2.7  is  analyzed  in  Appendix  B.  Equation 
(B-17)  from  this  Appendix  is  the  resultant  beam  element  matrix  equation 
neglecting  the  elastic  stiffness  contribution  of  the  actuator 


[K,] 

where  is  the  piezoceramic  actuator  coupling  coefficient  and  Vq  is  the 
applied  voltage  to  the  actuator  electrodes.  As  noted  in  eq.  (2.41),  the 
contribution  of  an  applied  DC  voltage  to  the  actuator  is  equivalent  to 
the  application  of  two  equal  an  opposite  couples  at  the  element  nodes. 

If  an  AC  voltage  is  applied  to  the  actuators  then  the  applied 
moments  would  be  dyr  mic  in  nature  and  the  response  of  the  FEM  model 
could  be  determined  given  the  mass  matrix  of  the  system. 

For  a  system  considering  only  the  mass  of  the  structure  and  not  the 
transducers,  the  resulting  equation  of  motion  could  be  written  as 
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I 
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In  reference  35,  the  same  element  eq.  (2.44)  is  obtained  by  using  the 
method  of  weighted  residuals.  When  the  ratio  of  the  thickness  of  the 
actuator  to  the  beam  is  much  greater  as  in  reference  35,  the 
contribution  to  stiffness  and  mass  of  the  transducer  must  also  be 
considered. 

The  structure  shown  in  figure  2.8  is  considered.  A  transducer  is 
bonded  to  the  upper  surface  as  a  sensor* similarly, the  transducer  bonded 
to  the  bottom  surface  is  a  actuator.  The  output  voltage  from  the  sensor 
is  increased  in  gain  by  a  constant  amount  and  supplied  to  the  actuator. 
Depending  on  either  the  orientation  of  the  polarity  vectors  of  the 
unimorphs  or  the  sign  of  the  feedback  voltage,  the  sign  of  the  couples 
at  the  element  nodes  is  determined.  With  no  other  type  of  conditioning 
of  the  sensor  signal  voltage  supplied  to  the  actuator  is  equal  to 


Vo(t)  =  VgCt)  =  GKg  I  ^ 


(2.46) 


where  represents  the  sensor  coupling  coefficient  of  eq.(2.20),  and 
X._j  and  X.  are  the  locations  of  the  element  node  points.  If  it  assumed 
that  the  applied  voltage  to  the  actuator  is  slowly  varying  such  that  the 
inertial  effects  of  the  beam  can  be  ignored.  The  resulting  beam  equation 
(Appendix  B)  for  a  collocated  sensor  actuator  is  written  as 


1  Sym.  e._ 


0 


(2.47) 


The  element  matrix  [C^q]  is  multiplied  by  the  velocity  of  the 
element  displacements  and  rotations.  As  such,  it  is  considered  as  an 
equivalent  viscous  damping  contribution  to  the  element.  This  matrix  can 
be  called  the  "electronic  damping"  matrix. 

A  four  element  approximation  of  a  cantilever  beam  with  various 
sensor  and  actuator  configurations  is  shown  in  figure  2.9.  By  adding 
the  additional  assumption  that  a  sensor  can  only  supply  a  feedback 
voltage  to  an  individual  actuator,  the  corresponding  global  system  [C^q] 
matrices  for  each  case  is  shown.  Case  B  of  figure  2.9  indicates  that  for 
non  collocated  sensors  and  actuators  that  the  [C^q]  is  no  longer 
symmetric. 

For  the  configurations  shown  in  figure  2.9  the  total  inertia  matrix 
for  the  complete  structure  is 

4 

CM]  =  [  [A^l^CM^lCA^]  .  (2.48) 

s=l 

where  [A]^  is  a  rectangular  matrix  that  is  used  to  assemble  the  element 
mass  matrices  into  a  mass  matrix  of  the  entire  structure. 

Similarly, 

4 

[K]  =  I  [A^l^CK^lCA^]  (2.49) 

s=l 

is  the  assembled  system  stiffness  matrix. 

Depending  on  the  global  location  of  the  sensor/actuator  active 
elements,  the  system  electronic  damping  matrix  is  also  assembled  as 


rjK  .'•ji  ’yji  -j*  ^  ~j^  '_«  ^  -*  ji  -_^  -p  -f  ^  -r  ^  ^  ---  j*  *  -.p  ^  ^  ^  -w  ^ 
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4 


[Cjd]  =  I  [B3]T[Cjp][B,] 
S  =  1 


(2.50) 


where  [B^]  is  a  rectangular  matrix  coordinating  the  location  of  the 
sensors  and  their  respective  actuators. 

The  vector  of  all  the  non  conservative  joint  forces  for  the  system 
with  electronic  damping  is  written  as 

4  4  • 

{fil=[[yV3)-V[8/[Cpp][B^]  h'j  .  (2.51) 

S=1  S=1  ’ 


The  equation  of  motion  for  the  complete  structure  is 
w. 


CM] 


V 


+  [K] 


6 


(Fi) 


(2.52) 


i, 


where  the  non  conservative  nodal  forces  are  defined  by  eq.(2.51). 
Combining  eqs.  (2.51)  and  (2.52)  and  excluding  all  viscous  type  forces 
from  the  non  conservative  force  vector  gives 


CM] 


1-1 

}  \ 

>  +  [C] 

r  • 

I  •  I 

'  hi 

}  I 

>+  [K]< 
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i  I 

hi 

(2.53) 


where  the  internal  damping  matrix  [C]  is  defined  as 

4 

CC]  =  y  CAJ^CC]CAJ  . 


(2.54) 


The  contribution  of  the  electronic  damping  matrix  in  eq.{2.53)  is 
directly  related  to  the  internal  system  viscous  damping  matrix  for  rate 
feedback  between  the  sensor  and  actuator  (no  filtering). 

It  will  be  shown  later  that  when  signal  conditioning  is  used  to 
selectively  increase  the  damping  for  a  particular  mode  that  for  other 
than  a  special  case,  the  feedback  will  add  to  both  the  stiffness  and 
damping  of  the  system. 

2.4.  Model  of  the  Feedback  Active  Control  Electronics 
A  specific  type  of  signal  conditioning  electronics  is  now  examined. 
Schematic  drawings  for  the  prototype  electronics  are  shown  in  Appendix  D. 
From  the  results  presented  in  appendix  D  the  final  model  of  the  active 
control  filter  unit  was  determined  to  be 


/q."  (VQ)  4  y  V3_ 


(2.55) 


where  is  the  filter  center  frequency 

Q  is  the  filter  bandwidth 

G  is  the  gain 

Vjj  is  the  voltage  to  the  actuator 

Vg  is  the  voltage  from  the  preamplifier. 


Solving  the  integral  of  Equation  (2.20)  for  a  particular  sensor 
element  yields 

Vs.{t)  =  Vo^^J  =  '^st^  *  ®i-l^  (2.57) 

Combining  equati ons(2 . 57)wi th (2 . 56)gi ves 

[I]{V  +  I^aJ  {Vp}  +  {Vjj}  =  (2.58) 

where 

[Tg]  is  the  sensor  location  matrix  containing  I's, 

-I's  and  O's  (rx  2n) 


For  the  structural  system  excited  by  only  actuators  the  matrix  equations 
are  written  as 

[M]  {q.}  +  CC]{q.}  +  [K]  {q.}  =  .  (2.59) 

The  actuator  force  vector  contains  only  nonzero  moment  terms  and  can  be 
expanded  to  yield 


(Fen)  =  [T„][KJ{V„} 


(2.60) 


where  [Tp]  is  the  actuator  location  matrix  matrix  containing  I's  ,  -I's 
and  O's  (2n  x  r)  and  [Rp]  is  a  diagonal  matrix  of  the  actuator  coupling 
coefficients. 

Combining  equations  2.56,  2.59  and  2.60  yields 

[I]{V  ^  +  [b]{VQ}  =  [K^][Tg]  {q,.}  (2.61) 

CM]  {q\}  +  [C]  {q.}  +  [K]  {q.}  =  [Tp][Rp]{Vp}  (2.62) 

Eqs.  (2.61 )  and  (2.62)  represent  the  coupled  pairs  of  differential 
equations  modelling  the  active  structure  conditioned  with  a  modal 
feedback  filter  system.  Previous  work  has  developed  a  similar  set  of 
equations  for  control  of  a  plane  grid  structure  using  accelerometers 
sensors  and  electrodynamic  shaker  actuators. 


There  has  been  a  significant  amount  of  reported  research  work  in 

the  field  of  the  identification  of  benign  structural  dynamic  systems  or 
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open  loop  systems  that  do  not  have  feedback  effects  However,  there 
is  no  reported  research  work  in  the  field  of  identification  of  an  active 
structural  dynamic  system  or  closed  loop  structural  dynamic  system  with 
feedback.  In  the  following  sections,  the  development  of  procedures  to 
identify  parameters  of  the  model  for  closed  loop  systems  and  the 
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3.1.  Equations  of  Identification 

A  system  shown  in  Figure  3.1  is  considered.  This  system  was  been 
approximated  by  n  finite  elements  each  with  a  unique  sensor/actuator 
pair  bonded  to  it  with  its  own  individual  signal  conditioning  system  or 
feedback  control  processor. 

Equations  (2.61)  and  (2.62)  in  Chapter  II  form  a  coupled  pair  of 
equations  describing  the  overall  feedback  control  system.  The  Fourier 
transform  of  Eq.  (2.61)  yields 

[-a,2[M]  +  ju)[C]  +  [K]]{q(a))}  =  {Fg(a>)}  +  {F^(oj)}  (3.1) 
where  as  before 


{q(cu)}  = 


w.  (u) 


{Fj.}  =  Active  actuator  control  moment  vector 


{Fg}  =  External  applied  forces  and  moments 


The  Fourier  transform  of  eq  (2.56)  yields 


where, 


[c.  -  a.o)  +  jwb^lVp. 


(u))  =  (w) 


(3.2) 


(3.3) 


a.  =  Q/Go)^,  b.  =  1/G,  c.  =  Qw^/G  (3.4) 

Taking  a  slightly  different  approach  than  Chapter  II,  the  control 
force  vector  {F^(a))}  when  expanded  has  the  following  form 

{F^(a))}  =  [0  (Mj  -  M2)  0  (M2  -  M3)  ...  (M^.J-  M^)  0  M^]^  (3.5) 

Solving  for  V^.  (u)  in  eq.  3.3  the  applied  moment  at  an  actuator  element 
node  is 


M^(u,) 


+  jwb. 


V.  (U)) 


(3.6) 


Using  eq.  (3.6)  and  the  fourier  transform  of  the  sensor  eq.  (2.20)  in 
Chapter  II  yields, 


M.(a))  =  Kp 


r - ^ - ^1[e(x.)-0(x.  )] 

i  [_  (c.  -  a.w  )  +  ju)b.  J  ^  ^  ^ 


(3.7) 


where, 


(3.8) 


By  using  a  notation 


hf  M  =  K  r - - 

i  i  |_  (c.  -  a.o)  )  +  jojb. 


(3.9) 


Therefore 


M^(u))  =  h^  (a,)[0(x.)  -  0(x..^)] 


(3.10) 


Equation  (3.5)  is  rewritten  as  follows 


{F^(aj)}  =  [0  0  -  [0  M^O  M3...M^0  0]^  (3.11) 


By  using  eq.  (3.9)  the  first  vector  component  is  equal  to 

ro  0  0  on 


T  0  h.  0  0 

[0  M,  0  M5...M j'  =  n 

^  "  0  0  0  0 


(3.12) 


■hi  0  hi 
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and  the  second  component  is  equal  to 


[0  0  M2..M^0  0]'  = 


0  0 

h-  0  -h, 
^2  2 
0  0  0 


Combining  eqs.  (3.12)  and  (3.13)  yields 


CFc(“)}  = 


where 


h-  +  h-  0 

Ti  Tg 


CH2(a.)]  = 


h^  +  h^ 


Combining  eqs.  (3.1)  and  (3.14)  further  yields 


{q^}  (3.13) 


By  denoting 


[H3((.)]  =  [[K]  -  a.2[M]  +  ju,[C]  -  [H2(a))]]'^ 


(3.14) 


(3.15) 


[[K]  -  +  jo,  [C]  -  [H2(a))]]{q(o))}  =  {Fg(u,)}  (3.16) 


(3.17) 


the  following  equation  is  obtained 


{q(ai)}  =  [H3(a,)]{Fg{a))}  (3.18) 

3.2.  Identification  of  Feedback  Filter  Electronics  Parameters 
Parameters  a.,  b.,  c.  of  eq.  (3,4)  can  be  calculated  from  a 
knowledge  of  the  circuits  of  the  signal  conditioning  system.  However, 
the  actual  values  used  for  active  control  were  determined  from 

experiments  which  measured  (t)  and  Vq  (t)  and  applied  a  least 

i  i 

square  identification  procedure  to  the  real  and  imaginary  spectral  data. 
The  admittance  function  of  equation  (3.4)  is 

Vs  ((o) 

^i^“^  "  V'p'^(Ji)  ■  *^i  (3.19) 

Expanding  eq.  (3.19)  into  real  and  imaginary  parts  yields 

Z,^^(a))  =  b.  (3.20) 

Zj^(a))  =  a. Cl)  -  c./u)  (3.21) 

The  real  and  imaginary  spectral  data  was  obtained  by  computing  the 
transfer  function  between  the  wideband  random  sensor  voltage  signal  and 
high  voltage  output  signal.  The  uncorrelated  noise  can  be  reduced  by 
computing  the  transfer  function  by  dividing  the  autospectrum  of  the 
input  voltage  by  the  cross  spectrum  between  the  input  and  output 
voltage. 

For  n  data  samples,  the  following  matrix  relation  can  be  written 
for  a  particular  feedback  system. 


0  -l/„„ 


(3.22 


where 


{y}  = 


a 
b  > 


c 

's_  ^ 


(3.23 


From  exact  values  of  a^,b.  and  c^,  equation  (3.19) is  exactly  satisfied 
For  other  values  of  a.,  and  c^  an  error  matrix  is  defined  as 
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{e}  =  {Z}  -  Cc]{y} 

and  a  square  measure  of  the  error  is 

1  _ 


(3.24) 


J  =  {e}'{e}  .  (3.25) 

Minimizing  J  with  respect  to  each  unknown  parameter  leads  to  the 
identification  of  the  parameters  a.,  b.  and  c^.  Each  sensor/actuator 
control  system  requires  n  independent  identification  computations. 

Equating  the  partial  derivative  of  eq.  (2.25)  with  respect  to  y  to 
zero  yields 


3^  =  0 

dy  3y 


(3.26) 

(3.27) 


Using  Eq.  (3.24)  and  eq.  (3.27)  solving  for  y  yields 
{y}  =  acficirhcf  {Z} 
Partitioning  the  matrices  of  eq.  (3.28)  yields 
{y}  =  CcJcj  +  cJ*C2]‘^{c{z^  +  cjzj} 
with  the  matrix  inverse  equal  to 


[cjcj  +  = 


ap  -  n 
0 


ap  -  n 
0 


(3.28) 

>.v, 

(3.29) 

i 

S:? 

**•  •  m 

(3.30) 

ap  -  n 


op  -  n 


Using  these  relationships  eq.  (3.29)  is  written  as 


{y}  =  -ib 
c 


J. 


n 


ap  -  n 
0 


n 


ap  -  n 


1/n 

0 


ap  -  n 

0 

CL 

ap  -  n^ 


d,}  (3.38) 


Solving  for  the  unknown  parameters  gives 


(d.p  +  d-h) 

a 

_  A 

2 

(3.39) 

ap  -  n 

b 

=  d2/n 

(3.40) 

nd,  +  d-a 

c  = 

i  0 

2 

(3.41) 

ap  -  n 

relations 

described  in  eq.  (3.4)  the  final 

results  are 

G  = 

1/b  gain 

(3.42) 

y  =  >/c/a  center  frequency 
Q  =  filter  bandwidth. 


A  computer  program  has  been  written  to  compute  the  feedback  system  gain, 
center  frequency  and  Q  from  experimental  admittance  data. 


The  procedure  for  identifying  the  feedback  filter  control 


coefficients  (a.,b.  and  c^)  is  as  follows.  A  wideband  1.0  volt  rms 
signal  was  sent  thru  a  lxlOE6  ohm  resistor  into  the  input  of  each 
control  unit.  This  produced  lxlOE-6  amp  wideband  input  into  the 
preamplifier.  With  the  gain  set  at  lxlOE6  the  output  voltage  was 
expected  to  be  l.OVrms.  The  transfer  function  was  taken  between  the 
input  voltage  and  the  output  voltage  with  the  control  units  set  for  the 
different  control  experiments.  The  frequency,  real  and  imaginary 
spectral  components  were  saved  and  processed  by  the  identification 
scheme  for  different  selected  bandwidths  about  the  filter  center 
frequencies.  As  a  comparison  to  the  identification  scheme,  the  center 
frequency  and  damping  (Q)  for  a  single  degree  of  freedom  model  curve  fit 
using  the  Genrad  complex  estimation  routine  51  was  used.  The 
identification  results  using  both  methods  are  listed  in  table  3.1. 

Figures  3.2  to  3.5  are  analytical  functions  using  equation  (3.7) 
and  the  identified  values  from  the  least  squares  method  over  the  actual 
experimental  wideband  transfer  function  data.  The  fit  for  the  first  mode 
is  better  than  the  second  mode.  The  high  gain  second  mode  fit  was  the 
poorest  (figure  3.5). 

3.3.  Baseline  Model  Development 

Equation  (3.28)  represents  the  model  of  the  structure  that  is 
controlled  by  using  piezoceramic  sensors  and  actuators  and  a  particular 
feedback  control  system.  The  approximately  known  feedback  control  system 
coefficients  (a.,  b.,  and  c.)  are  first  individually  identified  for  a 
I  particular  experimental  setting.  All  that  remains  to  be  identified  then 

I 
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gure  3-2  Comparison  of  Feedback  Filter  Analytical  Model  And 
Transfer  Function  Experimental  Data  -  Mode  1  Low  Gain 


i 


m 

I 

f* 

f 

f 

i 

I 


I  I  I  I 

•OVUI 


3  UJ 


a  >- 

<3 


%■ 


r 

r, 

? 


Figure  3-4  Comparison  of  Feedback  Filter  Analytical  Model 
And  Experimental  Transfer  Function  Data 
Mode  1  -  High  Gain 


is  the  dynamic  coupling  coefficients  defined  in  eqs.  (3.1)  and  (3.3)  and 
the  mass  stiffness  and  damping  matrices  of  the  baseline  non  controlled 
system. 

To  obtain  the  mass,  stiffness  and  damping  matrices  of  the  baseline 

system  an  identification  scheme  appropriate  for  structural  systems  with 
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non  proportional  damping  was  used.  A  brief  outline  of  the  procedure 
used  is  presented  in  Appendix  C. 

The  baseline  matrix  identification  system  requires  a  full  set  of 
eigenvalues  and  eigenvector  data  as  well  as  some  a  priori  mass, 
stiffness  or  damping  data.  From  this  information  the  identification 
scheme  is  used  to  identify  the  remaining  unknown  matrix  coefficients. 
Two  discrete  models  were  selected  for  the  cantilever  beam  shown  in 
figure  3.6.  Model  #1  consisted  of  five  translation  and  five  rotational 
degrees  of  freedom  and  model  #2  consisted  of  five  translational  degrees 
of  freedom.  The  two  models  were  identical  except  for  reduction  of  the 
rotational  degrees  of  freedom  in  model  #2.  Model  #1  represents  the 
coupling  of  the  transducer  in  the  structure  more  accurately  because  of 
the  retention  of  the  rotational  degrees  of  freedom  but  has  the 
limitation  of  the  difficulty  of  obtaining  rotational  eigenvector  data 
experimentally.  Model  #2  uses  Guyan  condensation®®  to  eliminate  the 
rotational  degrees  of  freedom  and  therefore  is  expected  to  increase  the 
overall  identification  error.  However  the  requirement  of  using  only 
translational  eigenvector  data  which  are  more  accurately  and  easily 
measured  justifies  the  investigation  of  using  such  a  model.  The  mass 
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stiffness  and  damping  matrices  for  model  #1  are  identified  using  the 

following  information: 

a)  selected  mass  elements  from  a  STRUDL®^  finite  element  model 

b)  experimental  eigenvalues  and  damping  for  modes  1  to  5 

(Appendix  0} 

c)  experimental  translational  eigenvector  data  for  modes  1  and  2 

d)  the  analytical  eigenvalue  data  for  modes  6  to  10  with  an 

estimated  .005  damping  value  for  these  modes 

e)  analytical  translational  eigenvector  data  for  modes  3  to  10 

f)  analytical  rotational  eigenvector  data  for  all  the  modes. 

For  model  #2  the  following  has  been  used: 

a)  selected  mass  elements  from  the  finite  element  model 

b)  experimental  eigenvalue  and  damping  for  mode  1  to  5  (Appendix 


c)  experimental  translational  eigenvector  data  for  mode  1  and  2 
(Appendix  D) 

d)  analytical  translational  eigenvector  data  for  mode  3  to  5 

The  identified  mass,  damping  and  stiffness  matrices  are  listed  in 
tables  3.2  and  3.3. 


To  verify  the  accuracy  of  the  identification  of  the  baseline 

matrices  a  computer  program  has  been  developed  to  compute  the  transfer 

function  matrix  of  either  the  5  or  10  of  models.  The  analytically 

generated  transfer  function  matrix  results  were  compared  to  baseline  non 

45 

controlled  experimental  results.  The  identification  scheme  can 
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Table  3.2  Identified  Baseline  Stiffness,  Damping,  and  Mass  Matrices 
For  Model  #1. 


(LOWER  TRIANGLE) 


STIFFNESS  MATRIX 


.1043E+05 

•.5710E*04 

■.873aE*04 

•.695aE*04 

•.2321E+04 

.2257E*04 

.3391E-04 

.a422E«-0S 

•.2ai0E*03 

.2013E*05 

•,5077E^03 

-.6046E<-04 

•.2772E^03 

'.433dE*04 


.1098E+05 
.a34SE>04 
.4027E*04 
.ia31E>04 
.  1040E*04 
.  190aE«-03 


,ia9aE*05 

.214aE*03 

,i4aaE*-05 

.6262E*04 

,4220E«-04 


.1398E«-05 

.7700E+04 

.1463E*-04 

-.1169E->-04 


,6aa8E*05 

.1427E-05 

,6031E*05 


. 1815E*05 
.2219E^05 


.6290E*03 

.1652E*05 

.22a6E*03 

.1193E*04 

.4a62£>03 

.2142E>04 


.2a7aE*04  .13iaE*04  -.2253E*05  .3a0aE»04 


.7976E*03 

.2223E*04 

.5310E*03 

.2674E*04 


.a871E*03 


.4635E*04  -.1253E+04 


.59a3E-03 

.4277E*04 


.2199E»04  -.5963E*03 


DAMPING  MATRIX. 


MASS  MATRIX 
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.27a7E-02 

-.3751E-03 

-.1301E-02 

-.4191E-03 

.ia22E-02 

.1216E-02 

.7339E-04 

.3377E-02 
- « 6964E-05 

.123aE-02 

ivlv 

-.2700E-04 

-.5177E-03 

-.2136E-02 

.5314E-03 

.5440E-02 

.1951E-03 

-.9349E-04 

-.3872E-03 

-.2444E-03 

-.8060E-04 

.8791E-03 

-.a227E-05 

.4a36E-02 

.4497E-06 

.5284E-03 

-.1334E-03 

-.3792E-02 

.4142E-03 

-.733aE-04 

.4771E-03 

-.4845E-04 

.1208E-02 

.1929E-03 

.2390E-03 

-.7764E-03 

-.2174E-03 

.5578E-04 

-.1041E-02 

.5208E-04 

.1281E-03 

-.1709E-04 

.7587E-03 

.5662E-04 

.50a4E-03 

-.3209E-04 

-.1492E-04 

-.ai23E-03 

.1352E-03 

.1886E-04 

-.1479E-03 

.4809E-03 

-.1527E-04 

.7362E-03 

.5153E-03 

.1254E-05 
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t 

.4497E-04 

.7122E-05 

.4962E-05 

.2045E-05 

.6332E-05 

-.2143E-05 

-.6569E-06 

.3163E-04 

.6267E-06 

.1840E-05 

'.Vs- 

.5562E-06 

-.2560E-06 

.57a3E-05 

-.1772E-05 

.2311E-04 
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V.Vi 

.1737E-07 

-.liaOE-06 

.1750E-05 

-.6884E-06 

.2230E-05 

.9130E-06 

.3080E-05 

.3a39E-04 

.5255E--06 

.1336E-05 

.3735E-07 

.4409E-05 

-.4061E-06 

-.a535E-06 

-.1018E-04 

-.8135E-07 

.5064E-05 

-.aaoiE-o7 

.1357E-06 

,70a7E-07 

-.1992E-06 

M 

-.ill4E-05 

.1195E-04 

-.1047E-06 

-.6727E-0S 

-.1127E-06 

.4287E-04 

-.1637E-06 

.1754E-05 

.9116E-07 

.2321E-06 

.6743E-03 

.89a5E-07 

-.3733E-05 

-.4543E-08 

.lllOE-04 

.2174E-07 

.5017E-05 

.63a6E-06 

-.2405E-0a 

y':-] 

TABLE  3.3  Identified  Baseline  Stiffness,  Damping  and  Maas 
Matrices  for  Model  #2 

STIFFNESS  MATRIX 

.6051E»04 

-.5227E*04  .7623E*04 

-.460dE<-03  -.52a6E«04  .146dE-^05 

.3063E+04  .80a3E*03  -.1157E+05  .1167E+05 

-.55a4E*03  .10aaE*02  .ia5aE+04  -.2124E+04  .4465E*03 

DAMPING  MATRIX 

.2715E-02 

-.1433E-02  .2241E-02 

-.59a7E-03  -.7677E-03  .376dE-02 

.9774E-03  -.4103E-03  -.3073E-02  .4042E-02 

-.1325E-03  .75a6E-04  .3727E-03  -.6770E-03  .3467E-03 

MASS  MATRIX 

.S127E-04 

-.14iaE-04  .3923E-04 

-.1045E-04  .2311E-05  .5055E-04 

.20iaE-04  -.5171E-05  -.4475E-04  .9477E-04 

-.3630E-05  .35a3E-05  -.2714E-05  .1040E-04  .2214E-04 


only  identify  the  mass,  damping  and  stiffness  matrices  to  within  a 
scalar.  Therefore  the  analytical  and  experimental  transfer  functions 
differed  slightly.  The  identified  matrices  were  scaled  to  match  the 
experiment  data  at  one  frequency  location  (380  Hz).  Comparison  of  the 
real  and  imaginary  parts  of  the  transfer  function  data  of  the 
experimental  and  baseline  models  for  modes  1  and  2  is  presented  in 
figures  3.7  to  3.10.  In  table  3.4  the  scaling  values  have  been  listed. 

Both  models  adequately  identified  the  first  two  modes  of  the 
structure.  The  greatest  variance  between  the  models  and  the 
experimental  baseline  results  occurred  in  the  real  part  of  the  first 
mode. 

3.4. Identification  of  Dynamic  Coupling  Coefficients 
The  coupling  coefficients  and  the  sensor/actuator  parameters  are 
contained  in  as  defined  by  eq.  (3.8  ),This  constant  when  expanded 
yields 


S,  "  ^s^d  ^*^3^  *^31  ‘^ll‘^F®3l'^d^^ 


(3.43) 


with  units  defined  as 

[C  0HM][NM/\/]  =  NM  sec  or  in  lb.  sec. 


Let  c  =  represent  the  dynamic  coupling  of  a  particular  pair  of 

collocated  sensor/actuator  to  the  elastic  beam. 

It  is  assumed  that  the  following  transfer  function  matrices  are 


available, 


The  Baseline  Model  #2  (5  D.O.F.)  And  Transfer 
imental  Data  -  Mode  1  -  No  Control 
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[Hgx(‘*»)]  experimental  data 
[H2{<*>)]  analytical  prediction. 

An  objective  function  is  defined  as  the  difference  between  the 
experimental  and  identified  analytical  prediction  and  is  written  as 
follows: 


(ij2  n  ‘^3  ^ 

J  =  I  Hgjj  (a,)  -  H3'^(a,)|2da>  +  |  M  (a,)|2du)  +  ... 

uij  k=l  ^  ^  0)2  k=l  ^  ^ 

(3.44) 

where  for  selected  frequency  bands  the  error  is  defined  as 
|e,|  = 

is  squared  and  summed  over  all  the  transfer  functions  available. 
Minimization  of  the  objective  function  is  then  done  with  respect  to  the 
unknown  parameters  and  the  parameters  identified  using  a  modified 
Marquardt  approach^^ 

For  a  structure  with  a  single  collocated  sensor/driver  pair  eq. 
(3.44)  reduces  to 
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Using  eq.  (3.42),  eq.  (3.43)  can  be  rewritten  as 
NPTS  N  FUNCTIONS 


(3.46) 
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where  W.  is  the  appropriate  weighting  matrix  factors  for  the 
integration  and  NPTS  are  the  number  of  points  used  for  the  integration, 
and  *  denotes  the  complex  conjugate. 

To  identify  the  unknown  parameters  (dynamic  coupling  coefficients) 
the  objective  functions  defined  in  eqs. (3.42)and(3.45)are  minimized  with 
respect  to  these  parameters. 

Briefly,  suppose  that  a  is  a  vector  which  makes  the  objective 
function  J(a)  a  minimum.  Let  a’  be  the  estimate  of  a  after  the  ith 
iteration.  As  shown  in  reference  78  ideally  the  i+1  iteration  should  be 

=  a’  -  p[R.]  q.  (3.47) 

where  p  is  the  step  size,  is  a  positive  definite  matrix  and  q^  is  the 
gradient  vector  of  the  objective  function.  When  Matrix  is  positive 
definite  the  direction  of  the  next  iteration  is  acceptable  (i.e., 
improves  the  minimization).  When  R^  is  the  identity  matrix  the  method 
is  referred  to  as  the  steepest  descent.  To  improve  the  iteration 
method  an  approximation  of  the  objective  function  at  the  minimum  can  be 
used.  An  approximation  of  the  objective  functions  behavior  at  a  minimum 
is  achieved  by  a  second  order  Taylor  series  expansion  of  the  objective 
function  about  the  i^^  iteration  value.  Equating  the  derivative  of  the 
second  order  expansion  with  respect  of  the  unknown  parameters  to  zero 


<.■  ^  k  ^  4  .■ 
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where  H.  is  the  Hessian  Matrix. 

As  mentioned  in  reference  76  various  methods  have  been  developed  to 
solve  the  iterative  scheme  defined  by  equation  (3. 48j.  The  particular 
method  used  for  the  identification  of  the  dynamic  coupling  coefficient 
was  the  modified  Marquardt  method. Two  particular  features  of  this 
method  are  that  an  approximate  Hessian  matrix  is  used  and  iterative  step 
size  is  modified  prior  to  each  iteration.  The  method  required  close 
estimates  of  the  dynamic  coupling  coefficients  for  successful 
convergence.  This  was  not  a  problem  since  the  range  of  the  dynamic 
coupling  coefficient  was  between  zero  (unbonded)  to  one  (perfectly 
bonded) . 

The  first  and  second  order  derivative  necessary  for  the  method  were 
computed  as  follows.  Taking  the  derivative  of  eq.  (3.46)  yields 

NPTS  N  ..k 

If  =  -2  I  I  <<  5r)  (3.A9) 

i=l  k=l 

where  from  Eq.  (3.17) 

[H3(a,)]  =  [[K]  -  co^CM]  +  ja,[c]  -  (3.50) 

where  [H2(a))]  is  the  contribution  of  the  signal  conditioning  system. 

From  the  following  identity 

[A]'^  [A]  =  [I] 


I 

I 

I 


(3.51) 


CA]-1  =  -  CA]-1[  [A]][A]'l 


Since  the  only  unknown  parameter  is  c  from  eqs.  (3.46)  and  (3.50) 


[Hj]  =  -  [Hjicfj  HjlCHj] 


I5  [H3]  =  [HjlC  HjJCHj] 


From  eq.  (3.9)  and  the  definition  of  any  particular  non  zero 
element  of  H2  is 

hi  (w)  -  - 2 -  (3. 

ij  (c.  -  a.co  )  +  ju)b. 


Therefore 


[H2]  -  ’  CH2(u.)] 


Eq.  (3.55)  becomes 


fc  ^^3] 


CHsJQ  H2^[H3] 


By  setting 


[H3]Qh2^  [H3] 


Eq.  (3.49)  yields 
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(3.60) 


When  more  than  one  collocated  sensor/driver  pair  is  used  the 
objective  function  becomes 


NPTS 
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NPTS  N 
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(3.61) 


i=l  h=l  i=l  k=l 

Then  for  each  unknown  coupling  parameter 


NPTS  N 


3J  _ 


3H 


k  NPTS  N 
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or  using  eqs.  (3.47)  to  (3.49) 


NPTS  N 


NPTS  N 


II  =  2  Re(e^*B5)  +  2  ^  W?'’  Re(e'f*  b'^)  +  (3.63) 


i=l  k=l 


i=l  k=l 


where 


[B.]  =  ^  [H3] 


(3.64) 


For  the  Hessian  matrix  (eq.  3.48)  the  approximation  of  the  second 
derivative  of  the  objective  function  J  is. 
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(3.66) 


It  can  be  shown  that  eq.  (3.66)  reduces  to 
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for  example,  for  two  sensor/dri ver  pairs 
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3.5.  Analysis  Procedure  and  Program  Verification  Results 
A  computer  program  (EDMFPL.FOR)  has  been  developed  combining  the 


analytical  model  of  the  system  and  the  multiple  sensor/driver 


identification  scheme.  The  inputs  to  the  program  includes  of  the 
particular  model's  mass,  stiffness  and  damping  matrices,  frequency 
bandwidths  and  transfer  function  element  locations  of  interest,  feedback 
control  filter  identified  parameters,  sensor/actuator  locations, 
weighting  matrix  for  the  numerical  integration  and  measured  transfer 
function  data  corresponding  to  the  locations  of  interest.  For  model  #2 
(5  dof),  the  FEM  condensation  matrix  has  also  been  provided.  In  figure 
3.11  data  flow  diagrams  for  the  various  programs  have  been  illustrated. 
Prior  to  identifying  actual  dynamic  coupling  coefficients  the  routine 
has  been  tested  by  using  analytical  generated  results  from  a  simulation 
program  (EDAT2.F0R). 

Estimated  dynamic  coupling  and  feedback  control  (a.,  b^  and  c.) 
filter  coefficients  for  mode  1  and  2  control  have  been  used  to  generated 
a  set  of  transfer  function  data.  These  tests  were  conducted  using  only 
a  single  transfer  function  which  was  the  tip  translational  element  ie, 
Hgg  (model  #1)  or  Hgg  (model  #2).  In  figures  3.12  the  identified 
results  around  the  first  and  second  mode  for  the  test  #1  are  presented. 
The  routine  usually  identified  the  unknown  parameters  exactly  within  10 
iterations. 

Two  set  of  experiments  have  been  conducted  on  a  cantilevered  test 
beam  as  shown  in  figure  3.6.  Actual  details  of  the  experiments  are 
discussed  in  the  Appendix  D.  In  both  experiments  the  control  of  the 
first  and  second  modes  of  the  beam  subjected  to  impact  loading  have  been 
discussed.  The  identified  feedback  control  filter  settings  for  both 
series  of  tests  was  previously  presented  in  table  3.1. 


In  table  3.5  the  identified  dynamic  coupling  coefficients  for  the 
different  tests  are  listed  and  reflex  the  changes  in  the  coupling  of  the 
collocated  sensor/actuator  pair.  For  a  perfect  bonding  condition  the 
coupling  would  have  a  maximum  value  of  1.00. 

Figures  3.13  to  3.14  represents  the  quality  of  the  matching  of 
the  model  with  the  identified  dynamic  coupling  factor  plotted  over  the 
actual  data  for  both  models. 

3.6.  Discussion  of  Identification  Results 
From  the  final  identified  dynamic  coupling  coefficients  presented 
in  table  3.5  various  observations  were  made.  The  collocated  sensor  and 
actuator  pair  with  its  filter  tuned  to  the  first  mode  was  twice  the  area 
in  size  as  the  pair  used  to  the  control  the  second  mode  and  the  dynamic 
coupling  coefficients  reflect  this.  Also,  the  dynamic  coupling 
coefficients  did  not  change  substantially  as  a  function  of  feedback 
gain.  The  lower  order  model  (5  dof.)  results  using  condensation  to 
remove  the  rotational  degrees  of  freedom  differed  in  magnitude  from  the 
10  dof.  model  but  approximately  followed  the  same  trends. 

The  identification  scheme^^’^®  using  the  objective  function  and  its 
derivatives  as  defined  in  eqs. (3.42),(3.6l)  and f3.65)was  used  to  identify 
the  dynamic  coupling  coefficients  separately  to  determine  if  their  was 
any  interaction  between  the  first  two  widely  spaced  modes.  As  expected, 
the  identification  results  did  not  change. 


Low  Gain 
High  Gain 


1827 

2205 


1026 

1030 


5604E*12 

3910E*12 


Model  #2 

Low  Gain 

.2348 

.0648 

.5170E*12 

High  Gain 

.3201 

.0712 

,4410E<-12 

Gain  ( .962)  Mode 


were  sensitive  to  errors  in  the  center  frequency  which  therefore  should 


be  a  key  parameter  to  be  varied  in  the  identification.  The  predicted 
dynamic  coupling  coefficients  were  changed  by  16  to  20%. 


To  obtain  an  additional  measure  of  the  accuracy  of  the 
identification  system,  the  square  root  of  the  error  squared  divided  by 
the  maximum  magnitude  within  a  selected  frequency  band  was  computed. 
The  error  is  defined  as  the  difference  between  the  experimental  and 
analytical  predicted  transfer  function  given  the  identified  dynamic 
coupling  coefficient  (eq.  3.42).  The  results  of  these  computations  are 
presented  in  table  3.6.  The  first  mode  low  gain  results  were  the  best 
and  the  second  mode  results  were  the  poorest. 

The  limitations  of  the  identification  procedure  were  primarily  due 


a)  Idenification  of  the  single  dynamic  coupling  coefficient  per 


mode. 

b)  Damping  estimates  and  mode  shape  rotational  coefficients  of 
the  baseline  model. 

c)  Single  degree  of  freedom  model  of  the  feedback  filter 
electronics. 


Table  3.6  Error  Computation  Resulta  Between  Measured  and 
Analytical  Prediction  Using  Identified  Dynamic 
Coupling  Coefficients 


Note:  Value  reflects  error  summed  over  the  bandwidth  divided 
by  the  maximum  magnitude  within  the  frequency  band. 
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CHAPTER  IV 


MODAL  AND  PARAMETER  IDENTIFICATION  OF  STRUCTURES  WITH  NON  COLLOCATED 

SENSORS  AND  ACTUATORS 


In  Chapter  III,  a  parameter  identification  scheme  has  been 
developed  for  an  active  structure  with  collocated  sensors  and 
actuators.  In  this  chapter,  the  case  of  non  collocated  sensors  and 
actuators  has  been  considered.  Two  active  structures  shown  in  figure 
4.1  is  considered.  Case  A  is  an  idealized  structure  with  the 
sensor/actuator  pairs  non  collocated  and  case  B  is  the  corresponding 
system  with  collocated  sensors  and  actuators.  It  is  assumed  that  the 
translational  degrees  of  freedom  have  been  removed  by  using  an 
appropriate  condensation  technique.  Furthermore,  it  is  assumed  that 
only  external  moments  are  applied.  For  the  purpose  of  illustration, 
coupled  differential  equations  for  the  structure  and  the  control 
electronics  are  written  for  a  system  with  four  degrees  of  freedom. 
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and  fourier  transforms  of  eqs.  (4.1)  and  (4.2)  yield 
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These  matrices  were  inverted  by  using  the  MACSYMA  progran>  and  the 
results  for  the  two  cases  are  shown  in  Tables  4.1  and  4.2.  The  upper 
left  quadrant  of  each  of  the  resultant  matrices  represents  the 
measurable  transfer  function  elements.  From  Table  4.1  it  is  observed 


that  the  the  upper  left  quadrant  is  no  longer  symmetric  and  reciprocity 

between  the  off  diagonal  elements  is  no  longer  maintained. 

If  the  system  contains  a  non  symmetrical  transfer  function  matrix 

an  equivalent  linear  model  of  the  system  using  mass  stiffness  and 

44 

damping  matrices  will  have  non  symmetric  matrices  .  Such  a  system  is 

non  self  adjoint  and  orthogonality  of  the  eigenvectors  is  not 

maintained.  The  eigenvectors  of  such  a  system  and  the  eigenvectors  of 

the  transpose  of  the  system  are  orthogonal  to  each  other  and  form  a 

biorthogonal  set.  In  addition  to  making  the  system  no  longer  self 

adjoint,  non  collocation  of  sensors  and  actuators  results  in  compounding 

the  control  stability  problems  associated  with  the  control  of  flexible 
47  73 

structures  ’  .  One  example  of  the  difficulties  in  designing  an  optimal 

control  regulator  using  a  quadratic  performance  index  with  non 

52 

collocated  sensors  and  drivers  has  been  discussed  by  Cannon  .  A 
control  law  for  a  flexible  robot  manipulator  arm  with  a  tip  position 
sensor  and  a  root  torquer  has  been  developed  with  the  objective  of 

reducing  the  final  position  placement  of  the  end  of  the  arm. 

Modal  identification  of  a  structural  dynamic  system  with  non 

symmetrical  matrices  from  the  measurement  of  time  histories  of  applied 
forces  and  accelerations  at  selected  locations  requires  a  method  to  post 
process  experimental  transfer  function  data  and  extract  the  biorthogonal 
eigenvector  parameters. 
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Such  a  method  has  been  developed  by  Potter  The  modal  analysis 

algorithms  currently  used  to  extract  modal  data  from  appropriate 

51 

transfer  functions  impose  the  assumption  of  symmetric  system  matrices 


A 


'.i-  'AK _ 
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Structures  are  rarely  modelled  as  non  self  adjoint  systems  .  A 

popular  approach  to  solving  these  type  of  problems  is  to  expand  the 

unsymmetrical  matrix  into  a  symmetric  and  antisymmetric  matrix.  The  use 
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of  a  non  self  adjoint  analysis  has  been  used  by  Hablani  This 
analysis  yields  uncoupled  pairs  of  first  order  equations  containing 
adjoint  eigenvectors. 

The  intent  of  the  efforts  here  are  to  use  the  developments  by 
53 

Potter  and  identify  the  modes  of  a  non  symmetric  active  system  with 
active  feedback  control  from  the  measurement  of  the  time  histories  of 
the  applied  input  and  system  response  output.  The  following  is  a 
slightly  simplified  theoretical  presentation  than  that 
presented  by  Potter  and  serves  to  enhance  the  overall  discussion  of  the 
final  results. 


4-1.  Theoretical  Development 
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The  equivalent  linear  system  is  given  by 
[MjgCqV  +  +  [KgjCq^}  =  {F.} 


(4.7) 


where  all  the  matrices  are  considered  general  real  matrices. 
By  denoting 


the  system  equations  are  written  as 

[AlCy^}  +  [B]{y.}  =  {z.}  (4.8) 

where  i—  -i  i—  — i 


CO]  CM]^ 

-CM]s  CO] 

CA]  = 

.  CB]  = 

(4.9) 

JMl,  CC]^_ 

CO]  CK]^ 

• 

By  assuming  eq. 

{y,-}  =  {Y}  e^^ 

Eq.  (4.8)  becomes 

CCB]  +  XCA]]  {Y}  =  {0}  (4.10) 

for  the  no  excitation  case. 

By  solving 

det  |B  +  AA|  =0 

2n  complex  roots  and  (eigenvalue)  a  corresponding  eigenvector  f’^are 
deterimined.  Therefore 

CCB]  +  X  rA]]{/}  =  {0} 


(4.11) 


Since  [B]  and  [A]  are  non  symmetric  it  is  shown  that  eq.  (4.11)  can 


also  be  written  as 


[[B]  +  ApCA]]^  {rP}  =  {0} 


(4.12) 


Using  eqs.  (4.11)  and  (4.12)  the  biorthogonal  relationships  are 


{rP}^[A]{'p''}  =  0 
{rPf  [B]{/}  =  0 


(4.13) 


From  the  following  external  forces 
{Z^(t)}  =  {Z.}  e^“^ 
and  a  solution  of  the  form 


{Z^(t)}  =  {Z.}  e*' 


(4.14) 


(4.15) 


and  the  additional  assumption  that  either  type  of  eigenvalue  can  provide 
a  solution  of  the  form 


I 

(Y,)  =  I  a, 


(4.16) 


it  can  be  shown  that 


{Y^(t)}  =  ^ 


f  {r''}^{Z.}{/} 


(jw  -  A^) 


(4.17) 


Cr]'[B]W  =  [A] 

where  [A]  is  the  diagonal  matrix  of  the  eigenvalues,  then 


±  jo)  /  1  - 


Cr]'[A]W  =  [I] 

By  substituting  for  Y.  and  Z.  in  eq.  (4.17) 


An  expansion  of  eq. 

(4.18)  yields 

2n 

r 

(jw  -  A  ) 

(,,>  =  1 

r=l 

r 

For  force  input  at 

k  and  displacement  response  at  i 

2n  „r  ,k 

«ik 

L  -  A 

r=l 

If  Ar,  =  rriji^  represents  the  residues  at  A  eq.  (4.20)  becomes 

(4.21) 


In  order  to  identify  the  eigenvectors  and  the  adjoint  eigenvectors, 
it  is  necessary  to  identify  and  A.,  from  the  measured  data  in  the 


or  when  separated  into  components 


'ZV'Zn 


M22/«2n 


For  n  modes  then 


r'  =  [a}i/a}„ 


A^  /A^ 
IZ'^ln 


(4.30) 


.  1] 


=  [A^  /A^ 


a2  /a2 


■  1] 


=  TA"  /a" 


^12'^^ln  •  •  • 


With  the  relationships  established  in  eqs,  (4.29)  and  (4.30),  the 
biorthogonal  eigenvector  coefficients  can  be  identified  by  normalizing 
the  row  and  columns  respectably. 

From  eqs.  (4.29)  and  (4.30)  it  can  also  be  shown  that  obtaining  a 
full  set  of  residues  experimentally  is  not  necessary.  Given  the  residues 
of  the  lower  triangle  of  the  transfer  function  matrix  for  a  given  mode, 
the  upper  half  values  can  be  computed  from 


A*"  =  (a’'  a’’  )/a’' 

''ik  ^''ii  \k^'''ki 


(4.31) 


(4.32) 


and  l,m  cannot  be  such  that  an  upper  half  value  is  chosen  prematurely. 


4.2  Experimental  Results 

A  cantilever  beam  shown  in  Figure  4.2  represents  a  structure  with  a 
single  pair  of  collocated  sensor  and  actuator  and  a  single  pair  of  non 
collocated  sensor  and  actuator.  The  feedback  signal  is  conditioned  to 
correspond  to  the  second  mode  frequency. 

The  transfer  function  data  taken  at  the  two  reciprocal  locations 
noted  on  Figure  4.2  with  feedback  control  are  shown  in  Figure  4.3.  From 
this  data  the  difference  between  the  transfer  function  is  very  slight. 
To  avoid  any  shift  in  second  mode  frequency  from  a  repositioning  of  the 
accelerometer,  small  masses  were  used.  Results  for  the  beam  subjected 
to  low  gain  control  using  a  collocated  sensor  and  actuator  are  shown  in 
Figure  4.4.  The  difference  between  the  two  transfer  functions  was  more 
than  expected  indicating  that  factors  such  as  transducer  misalignment 
and  incomplete  bonding  may  be  present.  Figure  4.5  is  the  transfer 
function  data  for  the  non  collocated  control  case  shown  in  figure  4.2 
at  the  same  gain  level.  The  greater  reduction  in  the  second  mode 
response  at  the  same  gain  is  because  the  actuator  is  now  twice  the  size. 
Nevertheless  the  difference  between  the  transfer  function  results  is 
greater. 

To  obtain  a  numerical  estimate  of  the  difference  between  reciprocal 
transfer  functions  (H..,  H..)  The  difference  between  the  real  and 
imaginary  components  was  defined  as 


e  =  |H.j(a,)  -  Hj.(a))| 

=  |x.j(a))  -  Xj.(a))  +  j{y.j^‘^^  ~ 


(4.33) 


(4.34) 
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where  x..  and  y. .  are  the  real  and  imgainary  part  of  H... 

I J  I J  <  J 

Therefore  the  normalized  difference  is 


Ur 


* 

e  e 


= 


(4.35) 


AC 


NC 


where  ACj^^.  is  the  sum  of  the  difference  results  for  the  no  control  case. 

The  results  are  listed  in  Table  4.3.  The  numerical  results  support 
the  above  observations. 

Shown  in  Figure  4.6  is  another  beam  with  the  feedback  filter 
centered  at  the  second  bending  mode  frequency  using  a  single  sensor  and 
single  actuator.  Also  shown  in  the  figure  are  the  transfer  functions 
measurement  locations  necessary  to  obtain  the  lower  triangle  of  the 
transfer  function  matrix.  Using  the  residues  computed  from  these 
transfer  functions  the  residues  for  the  remaining  elements  of  the  full 
transfer  function  matrix  have  been  computed  by  using  a  post  processor. 
The  post  processor  consists  of  expanding  eq.(4.3l)as  follows: 

For  each  mode  r, 


'ik 


.  Ki  ^ 


(4.36) 


where  U. .  is  the  real  part  and  V..  is  the  imaginary  part.  Expanding  eq. 

'  V  *  J 

(4.36)yields  for  a  given  mode 
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Figure  4-6  BlorChogonal  Mode  Identification  Experltnen 
Configuration  And  Test  Matrix 
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(4.37) 


(4.38) 


The  same  post  processor  computed  and  normalized  the  left  and  right 
eigenvectors  according  to  eqs.  (4.29)  and  (4.30)  with  and  without 
feedback.  The  coefficients  of  the  second  left  and  right  eigenvectors 
were  used  to  compute  a  variance  given  as 


I  I  I 

i=l  j=l 


’  ^ji  =  ^ 


(4.39) 


to  provide  a  relative  measure  of  the  degree  of  biorthogonality  for  the 
different  cases.  If  the  eigenvectors  are  orthogonal,  then  e  will  be 
very  small.  For  example,  a  four  degree  of  freedom  finite  element  model 
of  the  beam  with  no  control  gave  a  £  of  .206xl0E-3.  Experimental 
results  when  no  control  is  applied  yielded  a  value  of  .124.  This 
number  can  be  considered  an  equivalent  measure  of  the  noise  of  the 
experimental  measurement  as  well  as  the  residue  estimation  and 
eigenvector  coefficient  identification.  For  the  case  of  non  collocated 
control,  the  experimental  results  were 


CHAPTER  V 


GAIN  OPTIMIZATION  OF  THE  ACTIVE  STRUCTURE 


In  this  chapter,  an  application  of  optimal  control  theory  for 
active  damping  of  the  structural  dynamic  system  consisting  of  multiple 
collocated  piezoceramic  sensors  and  actuators  is  discussed.  The 
control  objective  is  to  bring  the  active  structure  form  a  perturbed 
initial  state  to  the  origin  in  an  optimal  fashion. 


5.1.  Background 

To  apply  the  concepts  of  optimal  control  to  the  active  structure 
a  brief  review  is  presented  of  a  method  of  obtaining  optimal  feedback 
gains  of  a  system  described  in  state  space  with  full  state  and  limited 
state  output  feedback  applied.  A  detail  description  of  this  method 
for  full  state  feedback  is  provided  for  a  in  a  variety  of  texts  on 
modern  control  theory  54,58,61,78 

Consider  a  constant  coefficient  system  defined  as 


x(t)  =  Ax(t)  +  Bu(t) 


(5.1) 
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where  is  the  initial  time.  To  optimize  the  control  of  such  a 


78 


system,  a  quadratic  performance  can  be  written  as 
^f 

J  =  I  j  (xV  +  U^Ru)  dt 
^0 

where  Q  is  an  arbitrary  positive  semidefinite  matrix  and  R  is  an 
arbitrary  positive  definite  matrix. 

Let  the  functional  L  be  defined  as 

^f 

L(x,x,u)  -  J  j  C(x^Qx  +  u^Ru)  +  A^(-x  +  Ax  +  Bu)]dt 


(5.3) 


(5.4) 


where  A  is  a  vector  of  lagrange  multipliers.  Equation  (5.4)  can  be 
rewritten  as 


^f 

L(x,x,u)  -  \  ^  (H(x,u)  -  A^x)dt 


(5.5) 


where 


H(x,u)  =  ^[x^Qx  +  u^Ru  +  A^{Ax  +  Bu)], 


Setting  the  first  variation  of  L  to  zero  yields 


6L  =  0  =  6x  +  —  6x  +  6u. 

3x  g.  au 


(5.6) 


(5.7) 


%  V  • 


r 


y.  ^  ^ 
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Using  Equations  (5.5),  (5.6)  and  (5.7)  gives 


6L  =  [  (|^6x  +  5u 

J^'8x  3u 


A  6x)dt  =  0. 


Integrating  the  second  part  of  Equation  (5.8) 


t,  t. 


f  f  T  T  1^ 

A  6x  =  A  6X  -  A  6x  dt 

trt  trt 


0  “0 


The  first  term  is  zero  if  '5x(tQ)  =  0  and  A(t^l 
Combining  Equations  (5.8)  and  (5.9)  yields 

^f 

6L  =  I  [(  15  +  a"^)  3A  +  (|5)6u]  dt  =  0. 
Equation  (5.10)  can  be  simplified  when  A(t)  ij 


15  ---A^  .A(V=0 


Then 


«L=  ||S«u<lt 


(5.8) 

by  parts  yields 

(5.9) 

is  chosen  to  be  zero. 

(5.10) 

;  chosen  such  that 

(5.11) 


For  arbitrary  6u  we  are  led  to  the  necessary  condition. 


Using  Equations  (5.6),  equation  (5.11)  and  (5.12)  become 

=  -Qx  -  A^x  (5.13) 

u  =  -  R'^B^X  (5.14) 

Equations  (5.1),  (5.2),  (5.13)  and  (5.14)  are  the  necessary 

conditions  for  optimality. 

Selecting  a  solution  in  the  form 

X(t  )  =  P(t  )  X  (t  ) .  (5.15) 

Then  differentiating  equation  (5.15)  and  substituting  the  results 
in  equation  (5.13)  along  with  Equation  (5.1)  and  (5.14)  yields  the 
Ricatti  equation 

^  =  -PA  +  PBR'^b''’p  -  Q  -  A^P,  P(t^)  =  0  (5.16) 

follows  from  the  choice  of  X(t^)  =  0. 

It  can  be  shown  that  P(t,t^),  t<t^,  is  stable  backwards  in  time. 
Therefore, 

lim  ^(t,t-)  =  |S(t)  =  0  (5.17) 

t-*®  ^ 

With  Equation  (5.17),  equation  (5.16)  becomes  the  algebraic  Ricatti 
equation  for  full  state  feedback. 

T  -IT 

PA  +  A  P  -  PBR  B  P  +  Q  =  0 


(5.18) 


The  corresponding  optimal  control  in  state  feedback  form  is 


u  =  -R’^B^Px  (5.19) 

where  P  is  the  positive  definite  solution  of  equation  (5.18). 
Uniqueness  is  guaranteed  if  the  pair  (A,B)  is  controllable,  are  closed 
loop  stability  is  guaranteed  if  the  pair  (A,C)  is  observable,  where  C 
is  any  matrix  such  that  Q  =  C^C. 


Optimal  Output  Feedback 


Now  consider  the  situation  where  perhaps  all  of  the  state  cannot 
be  measured  and  the  control  vector  is  proportional  to  the  output 
vector.  Such  a  system  can  be  defined  as 


X  =  Ax  +  Bu 

X  ^  r” 

(5.20) 

y  =  Cx 

y  ^  R 

(5.21) 

u  =  -Gy 

(5.22) 

where  matrices  A,B,C  and  G  are  constant  coefficient  matrices  and  the 
initial  state  is  unknown,  but  assumed  to  have  zero  mean  with  a 
variance 

ECXoxT  )  =  (5.23) 

Optimization  of  such  a  problem  has  been  formulated  by  Levine  and 

Athans^^  and  later  by  Mendel.®^  An  algorithm  has  also  been  developed 
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by  Moerder  and  Calise  to  solve  the  problem  with  an  objective 


Jj 


function  defined  as 


J  =  I  r (xV  +  u^Ru)  dt|  +  y(G) 
■'0  ' 


(5.24) 


where  y{G)  allows  for  a  penalty  to  be  placed  on  certain  element  of  the 


gain  matrix.  Using  Equations(5.12)and(5.22)  equation (5.20) becomes 


X  =  [A  -  BGC]x 


whose  solution  can  be  of  the  form 


(5.25) 


x(t)  =  [exp{A  -  BGC)t]xQ 


*(t,0)  =  exp(A  -  BGC)t 


equation  (5.24)  becomes 


(5.26) 


(5.27) 


where 


K  =  r «^(Q  +  C^G^RGOt  dt 

J  A 


(5.28) 


(5.29) 


when  Equations  (5,27)  and  (5.22)  are  substituted  into  Equation  (5.24). 
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Furthermore  equation  (5.28)  can  be  written  as 


J  =  tr(K  Xj  +  y(G) 


(5.30) 


where  tr  denotes  the  trace  of  the  matrix  and  the  matrix  K  satisfies. 


(A  -  BGC)\  +  K(A  -  BGC)  +  Q  -  cVrGC  =  0  (5.31) 

To  obtain  the  optimal  feedback  gain  matrix  a  lagrangian  is  defined  in 
reference  60  as 

L(G,K,X)  =  y(G)  +  trKX^  -  tr[(A  -  BGC)\  (5.32) 

+  K(A  -  BGC)  +  Q  -  cVrGC]X^. 

Taking  the  partial  derivative  of  the  lagrangian  with  respect  to  each 
variable  yields 


-B^KAC^  +  RGCAC^  +  |yq(G)  =  0 

(5.33) 

(A  -  BGC)A  +  A(A  -  BGC)^  +  =  0 

(5.34) 

(A  -  BGO^K  +  K(A  -  BGC)  +  Q  -  cVrGC  =  0 

(5.35) 

where  Yq(G)  is  the  gradient  of  y(G). 

This  result  obtained  by  Moerder  and  Calise  was  solved  using  a 
sequential  numerical  algorithm.  The  algorithm  has  the  added  feature 
that  certain  elements  in  the  gain  matrix  can  be  penalized  to  reduce 
their  influence  on  the  system.  This  feature  is  of  particular 
importance  in  the  active  control  of  structures  represented  by  finite 
elements.  In  the  next  section,  the  optimal  gains  are  determined  using 
this  algorithm  for  a  simple  structural  system  subjected  to  active  rate 
feedback  control  using  piezoceramic  sensors  and  actuators. 


no 


5.3  Rate  Feedback  Model 


The  active  structure  shown  in  figure  5.1  consists  of  two 
collocated  pairs  of  piezoceramic  sensors  and  actuators.  The  output 
from  each  sensor  was  only  multiplied  by  a  gain  and  applied  to  its 
corresponding  actuator. 

The  following  assumptions  are  made  for  this  analysis: 

a)  The  closed  loop  system  is  stabilizable  by  output  feedback. 

b)  The  piezoceramic  sensors  and  actuators  do  not  contribute  mass 


- 


or  stiffness  to  the  structure. 

c)  The  piezoceramic  actuators  have  no  internal  dynamic 
characteristics  (actuator  dynamics). 

The  governing  equations  of  the  active  system  were  transformed 
into  the  limited  state  feedback  equations  (5.20)  to  (5.22).  The 
structural  was  modelled  as 
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r 
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where  the  translation  d.o.f.'s  have  been  eliminated  by  a  condensation 
I  procedure  and  the  vector  {M.}  represents  only  active  feedback  control 

moments.  When  control  network  #1  and  #2  of  figure  5.1  condition  the 
sensor  signal  in  a  rate  feedback  manner  for  only  first  two  modes,  the 
output  voltage  to  the  actuator  is 
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located  Active  Rate  Control 


where  represents  the  sensor  piezoceramic  and  dynamic  coupling 


coefficients. 
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(5.38) 


Combining  Equations  (5.37)  and  (5.38)  yields 
=  [Kp]  [G]  [K^]  {0.} 


(5.39) 


Substituting  Equation  (5.39)  into  (5.36)  and  multiplying  both  sides  by 
the  inverse  of  the  mass  matrix  yields 


{0.}  =  -[M]'Hc]  {0.}  -  [M]'^[K]{0.}  +  [M]’^[Kp][G][Kp]  {©.} 


(5.40) 


Let  the  following  state  variables  be  defined  as 
T 


{X}  -  {X^...X2g}  “  {0^...0g  ,  0j^...0g} 

Equation  (5.40)  can  now  be  rewritten  in  state  space  form  as 


(5.41) 


{X,}  = 


[0]  [I] 

-CM]‘^[K]  -[M]"^[C] 


{x^}  +  [B]  {u} 


•  V  v^' 
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••/vS 
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T 


(5.42) 


The  mass,  stiffness  and  damping  coefficients  have  been  obtained 
from  the  identified  model  of  the  selected  beam  (Model  #1),  shown  in 
figure  5.1.  The  sensor/actuator  coefficients  were  from  previous 
presented  results.  Perfect  dynamic  coupling  has  been  assumed  and 
Table  5.1  lists  the  values  of  the  matrices  A,B  and  C. 

Optimization  Results 

The  procedure  for  using  the  algorithm  is  as  follows: 

a)  The  numerical  values  for  the  A,B,C,Q  and  R  matrices  provided 
as  input  data  to  the  program. 

b)  The  program  requests  initial  gain  matrix  and  penalty  values. 

c)  The  program  computes  the  open  loop  eigenvalues  which  provide  a 
check  on  the  input  data. 

d)  If  no  initial  gain  matrix  is  given  the  program  estimates  one. 


Table  5.1  Numerical  valuea  for  Rate  Control  Analyala 
Model  #1. 


0.0 

0.0 

0.0 

0.0 

0.0 

-2.334E*07 
5.653E*-07 
-3.016E«-07 
-4 .60dE'*’06 
6.401D»06 


3.465E»07  ~1.938E+07 
-1.204E»08  l.a80E<-08 


0.0 

0-0 

0.0 

0.0 

0.0 

2.734E*06 


0.0 

0.0 

0.0 

0.0 

0.0 

1.097E*06 


1.190E*08  -2.429E*08  1.656E»08 
-1.090E»07  1.425E+08  -1.799E»Oa 
-6.459E*06  -9.703E*07  1.772E*08 


4.236E^07  -4.915E+06 
1.656E*08  -1.290E+07 


1.799E»Oa  5.264E«-07 

1.772E*0a  -7.933E+07 


1 .0 
0.0 
0.0 
0.0 
0.0 

-2.442E*01 

2.422E*01 

5.013 

3.169 

-5.384 


0.0 

1.0 

0.0 

0.0 

0.0 

1 .405E*01 
-4.755E*01 


0.0 

0.0 

1.0 

0.0 

0.0 

1.241 


0.0 

0.0 

0.0 

1.0 

0.0 

-8.169 


1.919E*01  -1.528 


2.173E*01  -6.656E*01 


4.390 

7.757 


2.619E»01 

-2.545 


3.201E*01 

-5.146E*01 

3.750E«-01 


0.0 

0.0 

0.0 

0.0 

1.0 

-1.856 
2.264 
4.129 
1 .349E+01 

-4.332E*01 


0.0 

0.0 

0.0 

0.0 

0.0 

-1 .503E*03 
7.544E*03 
■1.475E*04 
1 .006E*04 
-a.576E*04 


0.0 

0.0 

0.0 

0.0 

0.0 

-3.937E*03 
2.0a7E*04 
-4.326E*04 
3. 130E»04 
-2.738E*04 


-4.074E«-02  4.074E*02 

0.0  0.0 


0.0 

-4.074E*02 


0,0 

4.074E’ 


e)  The  iteration  procedure  begins  and  the  attempts  to  obtain  an 


optima!  solution. 

f)  After  a  particular  solution  is  obtained,  interactive  analysis 
is  possible. 

Using  previously  identified  structural  and  sensor/actuator 
coefficients  gave  the  following  open  loop  natural  frequencies  and 
damping  values 


Rad/sec 

Damping 

367.9 

.0090 

2329.8 

.0039 

6534.5 

.0039 

12918.1 

.0023 

20752.1 

.0023 

The  R  matrix  was  an  identity  matrix  and  the  Q  matrix  was  defined  as 


Mi 


-•  V 


penalizing  all  the  angular  velocities  equally  and  giving  no  penalty  to 
the  angular  displacements.  With  a  null  initial  gain  matrix  the 
resultant  optimal  gain  matrix  is 


V 
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Gain  Matrix 

2.419D-03  -2.570D-04 

-2.290D-04  1.515D-03 

with  a  performance  index  of  .9000+06  and  natural  frequencies  and 
damping  values  of 


split  over  damped 

When  the  Q  matrix  was  an  identity  matrix  almost  identical  gain 
and  eigenvalue  results  were  obtained.  The  optimal  gain  distribution 
increased  the  damping  of  the  higher  modes  compared  to  the  lower  modes. 
This  was  expected  since  the  contribution  to  the  objective  function  by 
the  higher  modes  is  greater  because  for  velocity,  each  displacement 
variable  modal  contribution  is  multiplied  by  its  corresponding 
resonant  frequency.  When  the  off  diagonal  terms  of  the  gain  matrix  are 
penalized  in  y(G),  the  resulting  optimal  gain  obtained  was 
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Gain  Matrix 


2.500D-03  -2.530D-07 

-7.600D-07  1.740D-03 

with  a  slightly  higher  performance  index  of  .9013D+06  with  the 
following  resultant  natural  frequencies  and  damping  values. 


Rad/sec 

Damping 

367.9 

.0091 

2330.2 

.0047 

6535.6 

.0085 

12917.7 

.0268 

split 

over  damped 

These  results  did  not  differ  much  from  the  results  obtained  from  not 
penalizing  the  off  diagonal  gain  terms  in  the  gain  matrix. 

To  show  the  effect  of  the  optimal  gain  distribution,  the  state 
equations  Equations  (5.20-5.22)  were  solved  using  a  Runge  Kutta 
technique  for  the  open  and  closed  loop  control  configurations.  The 
system  angular  rate  state  variables  were  given  a  unit  angular  rate 
initial  condition  and  the  time  histories  of  system  response  obtained. 

The  output  response  y^lt)  in  volts  from  the  sensor  in  control 
loop  #1  of  the  baseline  no  control  configuration  for  50  msec  is 
shown  in  Figure  5.2.  The  response  signal  content  contains  the 
contributions  of  the  five  modes  of  the  structure  with  their  respective 


clamping  values.  When  active  rate  control  is  initiated  for  the  same 


set  of  initial  conditions,  the  control  voltage  to  the  actuator  u^^  is 
shown  in  Figure  5.3.  Comparison  of  Figures  5.3  to  5.4  demonstrates 
that  the  control  signal  is  the  exact  inverse  of  the  sensor  signal  and 
scaled  appropriately. 

Since  the  optimal  gain  distribution  shown  previously  increases 
the  damping  of  the  higher  modes  for  the  system,  the  predominant  signal 
response  due  to  the  first  two  modes  does  not  change.  In  Figure  5.5 
the  baseline  and  active  control  output  signals  from  control  loop  #1 
for  the  first  20  msec  are  superimposed.  The  observed  smoothing  of  the 
output  signal  during  active  control  is  due  to  the  increase  in  the 
damping  of  the  higher  modes. 

The  active  structure  (Model  #2)  shown  in  Figure  5.6  was  also 
analyzed  using  the  limited  state  feedback  approach.  As  before  the 
output  voltage  is 
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(5.45) 


and  the  control  moment  vector 


0 


(5.46) 


The  system  equation  retaining  only  the  rotational  degrees  of 


Figure  5-3  Octive  Control  Input  Ui 


freedom  is 


0.  =  -[M]'^[C]{0}  -  [M]‘1[K]{0.}  +  [M]'^[Kq][G][K^]{0.}  (5.47) 


and  the  state  variable  chosen  were 


{x}^  =  {0j,  02*  0^,  ©2}- 


(5.48) 


Equations  (5.45)  to  (5.48)  are  rewritten  in  a  state  space  form 
and  the  optimal  gain  matrix  with  off  diagonal  element  penalized  is 


6.957  D-03  1.697  D-04 
2.239 D-05  5.778  D-03 


when  the  Q  matrix  is  defined  as 


and  matrices  A,  B  and  C  are  given  in  Table  5.2. 


The  resultant  system  baseline  and  changes  in  modal  damping  values 


Rad/sec 


Damping 

Baseline 


Damping  at 
Optimal  Gain 


.0015 


Figures  5.7  through  5.10  indicate  the  first  mode  and  the  second 
mode  response  attenuation.  The  results  were  similar  to  model  #1. 


Figure 


Figure  5.8  Active  Control  Dlaplaceeent  Reaponae  1  aec 


CHAPTER  VI 


OPTIMAL  CONTROL  USING  DISTRIBUTED  SENSORS  AND  ACTUATORS 

In  a  recently  reported  work  Bailey  and  Hubbard  have  developed 
an  optimal  control  algorithm  for  a  system  consisting  of  a  cantilever 
beam  with  a  point  sensor  and  a  distributed  actuator.  The  actuator  was 
a  unimorph  piezoelectric  transducer  made  out  of  polyvinyl idene 
fluoride  polymer  (PVF)  and  covered  the  entire  upper  surface  of  the 
beam. 

The  objective  of  this  chapter  is  to  extend  the  work  of  Bailey  and 
15 

Hubbard  to  explore  control  strategies  by  using  distributed 
piezoceramic  unimorph  transducers  as  both  sensors  and  actuators. 

6.1.  Analysis 

The  structure  shown  in  Figure  6.1  ,  consists  of  a  cantilever 
beam  of  length  L,  with  constants  E,I  and  p.  Subdomains  of  the  beam 
are  covered  by  pairs  of  bonded  collocated  piezoceramic  unimorph 
transducers  polarized  in  the  z  direction.  The  transducer  on  the  upper 
surface  is  used  as  a  sensor.  When  conditioned  by  a  zero  impedance 
amplifier  it  was  shown  in  chapter  II,  eq.(2.22)  that  the  output 
voltage  for  a  constant  width  sensor  is 

X.  <  X  <  X.  (6.1) 


where  represents  the  sensor  coefficient.  The  transducer  on  the 
lower  surface  is  used  as  an  actuator.  As  shown  in  chapter  II,  eq. 
(2.43),  the  moment  induced  by  the  actuator  at  any  point  within  its 
domain  is 

2 

M(x,t)  =  El  ^  -  K|jVp(t)  X.  <  X  <  Xj  (6.2) 

3X 

where  as  before  represents  the  actuator  material  and  dynamic 

coupling  coefficients  and  is  the  applied  voltage  to  the 

transducer  electrodes.  When  the  sensor  output  is  amplified  and  applied 
to  the  actuator  without  any  filtering  (rate  feedback)  VQ(t)  is 

VQ(t)  =  G  Vg(t)  (6.3) 


where  G  represents  the  gain. 

From  eqs.  (6.1),  (6.2)  and  (6.3)  the  applied  control  moment 
becomes  for  x.<x<x. 


1 

M(x,t)  =  El  ^  -KpGK^  1^ 


X-Xj 


3X 


,J 


(6.4) 


x=x.- 


As  previously  discussed  in  Chapter  II,  the  distributed  piezoceramic 
actuators  can  be  modelled  as  external  point  couples  applied  at  the 
ends  of  the  domain  xi ,  xj.  Using  the  conventional  Bernoulli -Euler 
beam  analysis  and  assuming  that  the  transducers  add  neither  mass  or 
stiffness  to  the  beam  the  governing  differential  equation  for  the 
portion  of  the  beam  with  the  collocated  transducers  is 


El  ^ 


(x,t)  _ 


f  I” 

i,  3X 


x=x . 
J 


Iw 

3X 


=  0 


(6.5) 


for  xi<x<xj. 

As  noted  in  Chapter  II,  the  sensor  output  is  approximated  by  the 
difference  of  the  angular  velocities  at  the  transducer  ends,  and  is 
not  a  function  of  x.  Therefore  in  eq.  [6.5),  the  second  spatial 
derivative  makes  this  contribution  vanish  within  the  domain  xXx<Xj. 
Nevertheless,  the  control  moments  are  present  in  the  moment 
equilibrium  equations  at  the  transducers's  boundaries  (ie.,x  =  x.  and  x 

=  Xj). 

Let  the  domain  of  the  beam  be  subdivided  into  three  subdomains 
defined  as  follows. 


w^:  0  <  X  <  X. 

W2:  x.<  X  <  Xj 

W3:  Xj<  X  <  L 

The  governing  differential  equation  and  boundary  conditions  for  each 
subdomain  are  as  follows 

4  2 

3  W,  3^W, 

El  - j  +  pA  - i  =  0  0<x<x. 

3x  3t 

3w, 

w  .  n  ...  .  * 


(6.6) 


I 


j  D  S  I,  3X  |x=Xj  3X  |x=x.  J 

Then  the  expression  for  the  bending  moment  at  any  point  in  the  beam  is 

M{x,t)  =  El  ^  (6.10) 

3x 


with  the  following  boundary  conditions 
x  =  0  w  =  0 


|«  =  0 
3X 


at  X  =  L 


3^w 

3X^ 


3!w 

3x3 


=  0 


=  0 


and  the  strain  energy  of  any  point  of  the  beam  is 


I  “(^1 

dx 

(6.11) 

J  ^  3x^ 

The  resultant  effects  of  the  piezoceramic  actuators  are  treated  as 
externally  applied  moments. 


6.2.  Distributed  Parameter  Control 


To  derive  a  control  algorithm  for  the  system,  Lyapunov's  second 
method®^  was  used.  This  method  is  based  on  the  principle  that  if  a 


system  has  an  asymptotically  stable  equilibrium  state,  the  stored  energy 


within  system  when  displaced  decays  with  increasing  time  until  it 

assumes  a  minimum  value  .  A  direct  relationship  can  also  be  shown 

between  Lyapunov's  second  method  and  the  generalized  quadratic 
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performance  indices  used  in  designing  optimal  control  systems 

A  functional  is  chosen  which  is  the  sum  of  the  strain  and  kinetic 
energies  of  the  free  response  motion  of  the  beam^^. 

F=  |j[El[^]^pA  (6.12) 

0 

and  for  the  system  described  yields, 
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Using  this  functional,  the  objective  will  be  to  bring  the  system  to 

equilibrium  from  a  disturbed  state  in  a  minimum  amount  of  time. 

70 

Alternately,  as  discussed  by  Komkov  ,  the  objective  will  be  to 


obtain  a  control  actuating  force  which  in  this  case  is  a  control 
couple  pair  that  has  the  characteristic  of  decreasing  the  total  energy 
at  a  maximum  rate  for  each  time  interval  under  consideration.  Such  a 


at  ^  ^  ax 


Integrating  the  second  part  of  each  integral  twice  yields 
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The  integral  in  the  brackets  cancel  yielding 
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(6.16) 


3^W2(x.,t)  3^W2(x.,t)N  .  3^W3(Xj,);)  3^W3(Xj,t) 


Equation  (6.16)  can  be  further  simplified  using  the  moment 


ndary  conditions  at  x.  and  Xj  giving 
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Finally  using  the  following  boundary  conditions  at 


X  =  X. 


3Wj^  Bv^2 


3X  3x 


2  2 
3  Wj  3  W2 


3x3t  3x3t 


and  at 


x  =  x. 


3W2 

~3x  “  sF" 


2  2 
B  »2  3  Wj 


3x3t  3x3t 


Equation  (6.17)  reduces  to 


=  KoGK^[9(Xj,t)  -  8(x.,t)]‘ 


(6.18) 


where 


0  =  ^ 
3x3t 


Equation  6.18  represents  the  time  derivative  of  F  given  the 
feedback  voltage  as 


Vfjit)  =  GK5[9(Xj.t)  -  8(x..t)] 


(6.19) 
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where  the  difference  in  the  angular  velocity  between  the  two  nodes  can 
either  be  negative  or  positive.  For  negative  feedback  the  gain  value  G 
is  set  negative  and  eq.  (6.19)  is  now  equal  to 

Vjj(t)  =  -GK^C0(Xj..t)  -  e(x.,t)]  (6.20) 

To  maximize  -  (3F/3t)  the  control  moments  at  x.  and  x.  can  be  such 

J 

that  a  maximum  value  is  always  applied  opposite  to  the  motion.  This 
is  done  by  setting  Vp(t)  to 

Vp(t)  =  -GKgSgn[9(Xj,t)  -  0(x.,t)]  (6.21) 

where  the  sgn(  )  is  either  +1  or  -1  depending  on  the  sign  of  the 

difference  in  the  brackets.  Eq.  (6.21)  will  be  referred  to  as 

Lyapunov  control  and  when  eq.  (6.20)  is  used  it  will  be  referred  to  as 

rate  control . 

This  result  is  similar  to  the  one  obtain  in  reference  15,  however 
there  is  no  need  to  approximate  the  control  law  by  using  a  linear 
velocity  measurement  as  was  done  in  reference  15  since  it  is  directly 
measured  by  the  piezoceramic  sensor.  Secondly,  the  derivation  has  been 
extended  to  sensors  and  actuators  occupying  a  subdomain  of  the  beam 
rather  than  the  complete  region. 

6.3  Finite  Difference  Evaluation 

To  evaluate  the  performance  the  optimal  control  laws  discussed  in 

the  previous  section  an  explicit  numerical  finite  difference 
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technique  was  used.  Consider  the  structure  shown  in  Figure  6.2, 
consisting  of  a  distributed  piezoceramic  sensor  on  top  and  an  actuator 
on  the  bottom.  It  is  assumed  that  the  transducers  add  neither  mass 
nor  stiffness  to  the  structure. 


The  following  analysis  consists  of  two  parts:  open  loop  which 
considers  only  the  beam  with  no  active  control  and  closed  loop  control. 


Figure  6~2  Finite  Difference  Model 


In  the  open  loop  control  the  beam  is  given  an  initial  velocity 
distribution  and  only  internal  viscous  damping  of  the  beam  is 
modelled.  In  the  closed  loop  control,  active  control  is  applied  to 
the  structure  given  initial  velocity  conditions.  The  output  from  the 
sensor  is  amplified  and  used  to  excite  the  actuator.  Two  approaches 
for  modelling  closed  loop  control  have  been  examined. 

Open  Loop  Model  -  No  Active  Control 

The  dynamic  behavior  of  the  elastic  structure  without  active 
control  is  modelled  as 


=  V 

at  '' 


M  -  FT  O 
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1  a^M 


(6.22) 


(6.23) 


(6.24) 


where  w  is  the  displacement,  v  is  the  velocity,  M  the  moment  and  c  is 
the  viscous  damping  coefficient.  The  boundary  conditions  for  the 


cantilever  beam  are 


w{0,t)  =  0  ,  ^  ° 

ax'^ 


g  (0,t)  =  0.  ^  ° 

3X^ 


The  beam  has  been  discretized  as  shown  in  figure  (6.2)  for  the 
development  of  the  finite  difference  equations.  The  nodes  NX  +1  and 


NX  +2  are  additional  fictitious  nodes  necessary  for  computing  the 
shear  and  moment  boundary  conditions  at  the  end  of  the  beam. 

Finite  differencing  approximation  of  eqs.  (6,22)  to  (6.24)  are 
for  node  1=2  to  NX  considering  the  n+1  and  n  time  step 

(w^'^^d)  -  w"(I))/At  =  v^’d)  (6.25) 


v”'*'^d)  -  v”d)  ,  cv”d)_  1  Pm  ^d+i)  -  2M”d)  +  M^d-i)~[ 

At  m  ■  m  L  J 

(6.26) 

"""'(■I,-  '^"(■1  =  ^  -  2v""'(I)  *  v"*l(I-l)] 

(6.27) 

For  node  1  the  following  finite  difference  equation  are  applied 

w^’d)  =  0 

w'^'^^d)  =  0 
v^l)  =  0 
v"^l(l)  =  0 

(6.28) 

m"(1)  =  ^  w"(2) 
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To  compute  the  moments  from  eq.  (6.27)  requires  that  the  velocity 


of  the  additional  fictitious  points  be  known  from  the  boundary 
conditions  at  x=L  (node  NX).  The  finite  difference  equations  can  be 


Finally,  the  following  shear  boundary  condition  is  used  to 
compute  the  velocity  at  the  n  +  1  time  step  for  the  NX  node 

m"'^^(NX+1)  =  m"‘^^{NX-1)  (6.38) 

To  obtain  only  the  first  mode  response  of  the  beam  the  following 
velocity  initial  condition  was  chosen  for  the  structure. 

For  node  1 

v"(l)  =  0  (6.39) 

and  for  nodes  2  to  NX 

v'^(I)=AR[( .36705)[sin6jX-sinh3jX]-( .5)[cosPjX  -  coshp^x]]  (6.40) 
where  AR  is  the  desired  maximum  velocity  at  node  NX  (tip)  and  3^^=1.8751 
for  a  cantilever  beam  of  length  one  for  the  first  mode.  Given  this 
velocity  distribution  the  beam  response  contained  only  the  first  mode 
contribution. 

The  computation  sequence  was  as  follows: 

a)  Initial  conditions  eqs.  (6.39)  and  (6.40)  are  solved  for 
time  step  1  for  nodes  1=  1  to  NX  . 

b)  For  1=2, NX  velocity  form  eq.  (6.35)  and  the  displacement 
from  eq.  (6.37)  are  computed.  Note,  the  moment  for  1=  1  to 
NX  +1  is  zero  for  the  initial  time  step. 

c)  The  displacement  for  the  fictitious  nodes  and  are  computed 
from  eqs.  (6.32)  and  (6.33). 

d)  For  nodes  I  =  2, NX  the  moment  for  time  step  n  +  1  is  obtained 


e)  The  next  time  step  is  incremented  and  the  procedure  steps 


(b)  to  (e)  are  repeated. 

Closed  Loop  Analysis 

Two  finite  difference  approximations  have  been  examined  to  model 
the  closed  loop  feedback  response  of  the  system  under  free  vibration 
with  a  first  mode  velocity  initial  condition  (eq.  6.39  and  6.40).  For 
each  approach,  two  different  control  laws  eqs.  (6.19)  and  (6.21)  have 
been  evaluated.  The  sensor  output  was  described  previously  as 


v,(t) 


^  3v  J 


(6.41) 


where  V.(t)  is  the  output  voltage  of  the  sensor  and  K  is  the  sensor 
s  s 

coupling  coefficients. 

First  the  sensor  is  assumed  to  cover  a  portion  of  a  beam  and  is 
discretized  as  shown  in  figure  6.3.  The  sensor  output  is  proportional 
to  the  difference  in  the  angular  velocity  at  its  end  points  provided 
no  nodal  lines  exist  within  its  domain  of  the  modal  frequency  of 
interest.  The  voltage  output  of -the  sensor  can  be  approximated  in  a 
finite  difference  manner  as 


V3^^(q) 


■v""l(q+] 


■v"'^(p+J 


2ax 


(6.42) 


where  p  and  q  are  the  beginning  and  end  nodes  of  the  sensor.  When  the 
sensor  covers  the  entire  surface  of  the  cantilever  beam  the  following 


additional  condition  is  necessary  at  the  node  p  =  1 


(p+l)  -  v"*(p-l)  =  0 


(6.43) 


The  first  method  of  solution  (approach  #1)  for  the  active 
structure  covered  entirely  by  a  piezoceramic  sensor  and  actuator  is  as 
follows.  The  output  sensor  signal  is  conditioned  by  a  controller 
(rate  or  Lyapunov)  and  fed  back  into  the  system  to  the  beam  tip  as 

m'’‘'^(NX  t)  =  KjjVp(t)  (6.4 

In  terms  of  finite  difference  eq.  (5.44)  can  be  approximated  as 


w"  -^(NX+l)  -  2w"  •‘(NX)  +  w"  •‘(NX-1)  = 


(6.45) 


m"‘^^(NX+1)  -  m'’'^^(NX)  =  0 


(6.46) 


form  the  shear  condition  at  node  NX.  Combining  eqs.  (6.45)  and  (6.46) 
expended  in  term  of  displacements  and  solving  for  w  (NX+1)  and  w(NX+2) 
yields 


w"'^^(NX+l)  =  2w'’'^^(NX)  -  w"'^^(NX-l)  +  ^  KpVp(t) 


w"‘^^(NX+2)  =  3w"'^^(NX)  -  2w"'^^(NX-l)  +  KpVp(t) 


In  approach  #1  the  control  was  implemented  through  the 
displacement  and  velocity  computations  for  the  fictitious  nodes  NX+1 
and  NX  +2.  A  more  direct  method  (approach  #2)  is  possible  by  using 
eq.  (6.24)  at  node  NX  with  active  control.  From  eq.  (6.2),  eq.  (6.24) 


becomes 


M  (NX,t)  _  CT  l_v  (NX.t)  ,  .  ft 
at  "  3^2 


Approximated  eq.  (6.49)  in  finite  difference  yields 


m"'^^(NX)  =  m"(NX)  +  [v'^‘^^{NX+l)  -  2v'^'^^(NX) 

Ax^ 

rv^+l  .  yll, 

+  v'^'^^NX-l)]  +  AtKp  °  ° 


Since, 


eq. (6. 50} becomes 


V^Ct)  =  K^V^(t) 


m"'*’^(NX)  =  m"(NX)  +  [v”'^^(NX+l)  -  2v'’'^^(NX)  +  v"'^^(NX+l)] 

Ax'^ 

*  KqK,  [v"*i  -  vj: 

For  rate  control  for  approach  #1  and  #2, 

V„  =  -  K^VJt) 

The  Lyapunov  controller  for  approach  #1  was 


1.0, <  0 

Vp  =  Kj  sgn  (V^)  where  sgn  (V^)  =-  0.0, =  0  (6.54) 

[-1.0,Vg  >  0 

The  Lyapunov  controller  for  approach  #2  was  not  pursued  since  it  is  a 
discontinuous  function.  Because  of  this  limitation,  approach  #1  was 
used  for  comparison  between  the  performance  of  the  rate  and  Lyapunov 
control  laws. 


The  computational  sequence  was  as  previously  except  that 
approach  #1,  step  C  now  uses  eqs.  (6.47)  and  (6.48)  and  either  eqs. 
(6.53)  or  (6.54)  for  control.  For  approach  #2  step  d,  for  I  =  NX  eq. 
(6.52)  is  used  and  eq.  (6.53)  for  control. 

Stability  Analysis 

To  compute  the  desired  response  of  the  finite  difference  model  of 

the  active  structure,  a  time  step  must  be  selected  which  insures 
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numerical  stability.  To  obtain  a  proper  time  step,  a  Von  Neumann 
stability  analysis  was  used  for  the  open  loop  control  which  in  turn 
served  as  an  upper  limit  to  empirically  find  the  proper  time  step  for 
the  closed  loop  control  model. 

The  two  coupled  finite  difference  eqs.  (6.26)  and  (6.27)  can  be 
written  as 


.,n+l  .  .,n  At  /^n  -^n  .  ^.n  , 


(6.55) 


un+1  _  „n  AtEI  ,..n+l 

^+1  -  ^  i"  rr  (''i 


AX 


(6.56) 


where  the  subscript  (i)  represents  the  spatial  node  and  the 
superscript  (n)  represents  the  time  increment.  The  fourier  component 
at  each  node  is  defined  as 


V 


=  v"ei«' 


,n,jei 


(6.57) 


mJ+I  =  M^+l  gj9(i±l) 


Using  the  relations  given  in  eq.  (6.57)  ,  eqs.  (6.55)  and  (6.56) 
become 


,n+l  =  ^  „n  ^  g-je^ 

Max'^ 


„n+l  =  .  At  ^n+l^^je  .  2  g-je^ 

LyT 


(6.58) 


(6.59) 


further  by  defining 


p  =  4  sin^  0/2 


Equations  (6.58)  and  (6.59)  can  be  further  simplified  and  written 
in  matrix  form  as 


t-  r  it-'i  ^  ii 
^  I  ..2J  " 


(6.60) 


The  eigenvalues  of  the  matrix  given  in  eq.  (6.60)  can  be  obtained 
by  solving  the  following  determinant 


1  -  A 


^  p 
2  ^ 
MAx'^ 


fM-f 

UjbL] 

M  J 

=  0  (6.61) 


>> 

X-*- 


I 


The  resulting  quadratic  equation  becomes 


A  -  A(2  -  A*)  +  1  =  0 


(6.62) 


Mhere 


2  ...2 


(M  T 


(6.63) 


The  roots  of  equation  (6.62)  are 


1  / - 2 - 

2  *  2  /  {Z-k*r  -  4 


(6.64) 


In  order  for  the  numerical  computations  to  be  stable  the  errors 
generated  must  not  grow.  It  can  be  shown  that  this  condition  is 
satisfied  when  if 


A  =  1  when  (2-A*)"^  =  4 
A  <  1  when  (2-A*)^  <  4 
A  >  1  when  (2-A*)^  >  4 


(6.65) 

(6.66) 

(6.67) 


Using  eqs.  (6.63)  and  (6.65)  or  (6.66)  the  time  increment  can  be 
directly  computed  given  the  spatial  mesh  size  from 


2  2 
W 

At  <  - - 

-  2a,. 


(6.68) 


where 


CO.  is  the  highest  natural  frequency  of  interest, 


p.  is  the  coefficient  for  the  highest  eigenvalue  which  is 
dependent  on  the  boundary  condition, 


P=1  for  the  maxitnum  value. 


Closed  loop  stability  was  evaluated  empirically  and  was  an  order 
of  magnitude  smaller  than  the  open  loop  time  step. 

6.4  Discussion  of  Results 

Optimal  Control 

A  seven  node  model  of  the  beam  shown  in  figure  6.2  was  used  for 
evaluation  of  the  different  models  and  control  laws. 

To  check  on  the  accuracy  of  the  active  structure  model  for  the 
first  mode  response  the  beam  was  given  the  first  mode  profile  initial 
condition  velocity  with  a  maximum  tip  velocity  of  1.0  m/sec  and  allowed 
to  oscillate  with  no  control.  Figure  6.4  is  a  plot  of  the  tip 
displacement  versus  time  for  this  test.  The  first  natural  frequency  of 
the  model  was  approximately  5.2  Hz  which  is  the  theoretical  natural 
frequency  for  the  structure. 

Figure  6.5  is  the  tip  displacement  response  given  the  same 
initial  condition  of  the  beam  with  internal  viscous  value  (c)  of  .002. 
Figures  6.6a  and  6.6b  are  the  feedback  actuator  signal  and  the  tip 
displacement  respectively  using  approach  #1  for  a  feedback  gain  of 
.036.  Figure  6.7  is  the  tip  displacement  using  approach  #2  for  the 
same  feedback  gain  and  initial  condition.  Figures  6.6b  and  6-7  are  in 
close  agreement. 

Figures  6.8a  and  6.8b  are  the  same  results  employing  the  Lyapunov 
controller.  The  controller  required  a  scaling  value  of  the  sensor 
maximum  output  for  the  given  initial  velocity  condition.  The  Lyapunov 
controller  was  disengaged  when  the  sensor  output  was  less  than  node 
separation  distance  times  the  maximum  initial  tip  velocity.  Comparison 
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of  figures  6.6b  and  6.8b  demonstrates  the  enhanced  control  performance 
of  the  Lyapunov  controller  using  piezoceramic  sensors  and  actuators. 
Sensor  Size  Evaluation 

Either  type  of  feedback  control  is  directly  proportional  to  the 
input  sensor  signal.  As  the  sensor  size  is  decreased,  the  sensor  output 
is  reduced  since  the  difference  between  the  angular  rate  at  the  two  end 
nodes  decrease  as  the  nodes  approach  each  other.  This  can  be 
demonstrated  using  approach  #2  to  evaluate  various  sensor  sizes. 

Figure  6.9  is  the  tip  displacement  of  the  rate  feedback  control 
using  a  .2  gain  level  with  a  sensor  covering  the  entire  upper  surface. 
Using  the  same  feedback  gain  level,  the  sensor  was  decrease  72%  in  size 
which  resulted  in  a  decrease  in  control  performance  as  shown  in  figure 
6.10.  Figure  6.11  demonstrates  that  control  can  performance  can  be 
reclaimed  by  increasing  the  feedback  gain  by  a  factor  of  18.5. 

As  mentioned  in  Chapter  II,  the  sensor  signal  will  degrade  for 

frequencies  at  which  there  exist  standing  nodal  line  within  its 

subdomain.  As  the  sensor  size  is  decreased,  the  bandwidth  of  the 

measured  input  to  the  control  system  is  increased  but  at  a  cost  of  a 

decrease  in  the  signal  to  noise  of  the  input  current  signal.  This  is 

usually  not  a  problem,  though,  since  with  the  proper  signal 

conditioning  electronics  the  sensitivity  of  the  piezoceramic  dynamic 
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strain  gages  is  very  high  . 

Actuator  Size  Evaluation 

A  similar  difficulty  arises  in  the  actuator.  As  noted  before,  the 
piezoceramic  actuator  is  modelled  as  a  system  with  dynamic  opposing 
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couples  at  its  boundaries.  Should  a  transducer  attempt  to 
control  a  mode  such  that  it  results  in  at  least  one  node  line  existing 
within  its  subdomain,  a  degradation  of  performance  would  result 
because  of  the  phase  difference  due  to  the  spatial  distance  between 
the  two  dynamic  couples.  To  eliminate  this  problem  requires  proper 
conditioning  of  the  control  voltage  to  remove  the  signal  contributions 
above  the  control  bandwidth  of  the  system. 

Incorporation  of  the  feedback  filter  electronic  differential 
equation  discussed  in  Chapter  II  into  the  finite  difference  scheme  was 
attempted.  Stability  difficulties  were  encountered  which  were 
primarily  due  to  the  second  order  differentiation  of  the  control  signal 
Finite  difference  evaluations  were  conducted  on  a  model  that  contained 
an  actuator  which  occupied  a  subdomain  of  the  beam.  Stability 
difficulties  were  again  encountered.  The  primary  suspected  problem  was 
the  inability  to  adequately  model  the  point  nodal  moments  in  the 


domain  of  the  beam. 


CHAPTER  VII 


ACTIVE  CONTROL  EXPERIMENTS 


Throughout  out  this  research  effort,  a  variety  of  active  control 
or  "electronic  clamping"  experiments  were  undertaken.  The  experiments 
consisted  primarily  of  investigating  the  behavior  and  performance  of 
the  individual  unimorph  transducers  or  analysis  of  the  behavior  of  an 
active  structural  system.  The  experiments  were  inherently  flexible 
due  to  the  adjustability  of  the  feedback  control  electronics,  to 
explore  both  rate  and  modal  feedback  system  stability  and  performance 
issues. 


7.1.  Experimental  Configuration 

A  typical  test  configuration  is  shown  in  figure  7.1  A  test 
structure  was  chosen  and  unimorph  transducers  bonded  to  select 
locations.  In  all  cases  a  cantilever  beam  was  the  test  structure. 
Depending  on  the  desired  test,  transducers  served  as  either  sensors  or 
actuators.  The  output  of  a  sensor  was  conditioned  by  the  feedback 
electronics  which  in  turn  supplied  a  high  voltage  control  signal  to  an 
actuator.'  For  the  tests  conducted  in  this  research  one  sensor 


controlled  an  individual  actuator.  It  is  possible  to  control  more 
than  one  actuator  with  a  single  sensor  but  unfortunately  the  prototype 
electronics  were  electrically  unstable  in  this  type  of  configuration. 


Typical 


The  feedback  electronics  were  specifically  designed  to  be  used 
with  piezoceramic  transducers.  Each  unit  basically  consisted  of  a 
preamplifier,  filter  and  high  power  output  section.  The  preamplifier 
converted  the  current  output  of  the  sensor  to  a  voltage  signal  which 
was  conditioned  by  a  bandpass  filter  whose  center  frequency  and  Q  were 
adjustable.  The  conditioned  signal's  phase  could  also  be  adjusted. 
Finally  the  low  voltage  signal  was  amplified  and  applied  to  an 
actuator.  Details  of  the  different  experimental  configurations  and 
the  feedback  control  electronics  is  provided  for  in  appendix  D. 


Mo  da 1  Control  Experiments 


The  experimental  setup  provided  an  ideal  platform  to  examine  some 

of  the  problems  associated  with  the  application  of  modern  control 

theory  to  flexible  structures.  A  major  difficulty  when  this  is 

attempted  is  that  flexible  structures  are  distributed  parameter 

38 

systems.  According  to  Balas  the  fundamental  problem  of  feedback 

control  of  flexible  systems  is  meeting  the  requirement  of  precise 

control  of  a  large  dimensional  system  with  a  much  smaller 

dimensional  controller.  The  approach  used  by  most  investigators 
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reduces  to  the  problem  of  modal  control  Modal  control  as  defined 
40 

by  Simon  as  a  method  of  control  which  changes  the  eigenvalues  and 
eigenvectors  of  the  system  matrix  to  achieve  desired  control 
objectives. 


When  only  a  limited  number  of  modes  of  the  infinite  number  of 
modes  of  the  actual  structure  are  controlled  observation  spillover 


i 


usually  occurs.  Observer  spillover 


results  when  a  finite 


dimensional  controller  attempts  to  control  an  infinite  dimension 
system  such  as  an  active  structure.  At  some  stage  the  sensors  of  the 
system  will  sense  the  modes  of  the  system  that  are  not  accounted  in 
the  analysis  (residual).  This  will  lead  to  system  instabilities 
observed  when  the  narrow  bandpass  filters  are  bypassed.  Control 
spillover  occurs  when  in  the  process  of  controlling  the  desired  modes 

the  uncontrolled  modes  are  excited. 
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According  to  Balas  control  spillover  only  causes  a  degradation 

of  system  control  performance  while  the “observer  spillover  can  make 
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the  system  go  unstable.  Meirovitch  has  shown  that  even  a  small 

amount  of  structural  damping  in  the  uncontrolled  modes  has  the 

potential  to  eliminate  the  instability  caused  by  observer  spillover. 

The  use  of  the  active  control  feedback  with  arrow  bandpass 

filters  reduces  the  effect  of  the  residual  modes  on  the  control 

system.  Since  piezoceramic  unimorph  sensors  are  sensitive  to  very 
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small  dynamic  strains  they  tend  to  be  overly  affected  by  modes 
outside  of  the  control  bandwidth  and  thus  required  the  bandpass 
filters. 

Examples  of  theoretically  and  experimental  results  in  modal  space 
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control  are  presented  by  Meirovitch  ,  Hallauer  ’  just  to  name  a 
few.  In  these  works,  modal  space  control  has  been  accomplished  by 
using  fewer  sensors  and  actuators  than  control  modes  and  employing 
narrow  band  filters  similar  to  the  ones  used  in  these  experiments. 

To  initiate  rate  control  the  bandpass  filter  could  be  bypassed. 
When  this  was  done  experimentally,  as  expected,  the  control  feedback 


system  went  unstable.  Therefore,  all  the  experiments  used  the 
bandpass  filters.  Since  the  filters  only  had  a  Q  of  20  instabilities 
would  occur  if  the  current  to  voltage  gain  was  increased  passed  10^. 

To  measure  and  reduce  the  dynamic  response  data  a  16  channel 
Computer  Aided  Test  System  (Genrad  Model  2515)  was  used.  This  system 
had  spectral,  frequency  response  and  modal  analysis  capabilities. 

7.3.  Test  Procedures 

The  primary  form  of  testing  used  to  evaluate  the  modal  control 
performance  was  impact  response  transfer  function  analysis.  The  test 
procedure  followed  consisted  primarily  of  measuring  the  structure's 
response  from  a  stationary  accelerometer  while  impacting  the  structure 
with  a  force  gage  hammer  at  selected  locations. 

From  this  data  changes  in  the  system's  eigenvalues  and 
eigenvectors  could  be  extracted  using  the  various  methods  provided  for 
by  the  SDRC  Modal  Plus  software^^.  Prior  to  evaluating  the  system's 
response  due  to  impact  excitation,  the  center  frequency  of  the  active 
control  electronic  had  to  be  tuned.  Setting  the  center  frequency  at 
measured  frequencies  of  the  structure  did  not  provide  maximum  modal 
control  and  an  iterative  procedure  was  required  to  optimize  the  tuning 
of  the  filters.  When  the  gain  was  changed  the  filters  required 
retuning  to  maintain  maximum  performance. 


Figure  7.2  is  typical  of  the  data  obtained  during  the 
experimental  investigation.  Shown  in  the  figure  is  multimode  control 
of  the  first  and  second  bending  modes  of  a  cantilever  beam  at  two 
different  gain  levels  using  two  collocated  pairs  of  piezoceramic 
sensor  and  actuators.  The  low  gain  results  are  at  the  lowest  gain 
settings  of  the  feedback  electronics.  Both  high  Q  resonant  peaks  are 
reduced  with  the  relative  change  greatpr  in  the  higher  mode.  As  the 
gain  is  increased  to  just  below  the  point  at  which  the  system  goes 
unstable  the  peaks  are  further  reduced.  At  this  level  the  first  mode 
is  also  beginning  to  split.  Appendix  D  provides  the  actual  numerical 
data  for  these  particular  results. 

The  results  from  initial  tests  conducted  on  a  steel  beam  (Beam  #1 
appendix  D)  were  used  primarily  in  a  preliminary  study  of  identifying 
the  changes  in  the  damping  matrix  of  a  discrete  model  of  the  test 
cantilever  beam.  These  results  are  detailed  in  reference  77  and  the 
change  of  selected  damping  matrix  elements  were  determined  to  be 
proportional  to  the  gain  for  both  single  and  multimode  control. 

The  next  series  of  tests  examined  collocated  and  non  collocated 
control  of  a  aluminum  beam  (Beam  #2  appendix  D)  and  these  results  were 
discussed  in  chapter  IV.  The  last  series  of  tests  were  conducted  on 
the  smallest  test  structure  (Beam  #3,  appendix  D)  to  identify  the 
dynamic  coupling  coefficients  of  the  unimorph  transducers.  The 
results  of  these  efforts  is  presented  in  chapter  III. 


Rasellne.  Low  And  HIeH  Gain  Comparisons 


CHAPTER  VIII. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  potential  of  using  piezoceramic  transducers  for  active 
control  of  structures  has  been  demonstrated  in  the  past  by  a  variety 
of  experiments.  Based  on  the  results  of  these  efforts  the  objective 
of  this  research  has  been  to  understand  how  the  use  of  piezoceramic 
sensors  and  actuators  can  actively  control  a  flexible  structure  by 
developing  suitable  models  and  conducting  experimental  investigations. 

A  finite  element  model  of  the  active  structure  was  developed  from 
the  basic  piezoceramic  constitutive  equations  for  a  beam  element  with 
a  piezoceramic  unimorph  sensor  and  actuator  bonded  to  it.  From  this 
model  the  current  output  of  the  piezoceramic  sensor  was  determined  to 
be  proportional  to  the  integral  of  the  strain  rate  over  the  element 
length.  The  effect  of  the  actuator  was  approximated  by  a  pair  of 
opposing  couples  applied  at  the  element  modal  points  which  were 
proportional  to  the  applied  actuator  voltage.  Furthermore,  it  was  also 
established  that  no  collocated  sensing  and  actuating  transducers 
resulted  in  a  non  self  adjoint  system. 

Using  a  frequency  formulation  of  the  coupled  differential 
equations  of  the  active  structure,  experimental  data  and  least  squares 
type  of  identification  routines,  the  identification  of  the  dynamic 
coupling  coefficients  for  a  cantilever  beam  subjected  to  multimode 


control  was  accomplished.  Improvement  in  the  values  of  the  identified 
coefficients  is  possible  based  on  the  accuracy  of  the  baseline  and 
active  control  feedback  electronics  models. 

A  study  was  undertaken  to  examine  the  changes  in  the  active 
structure's  modes  when  controlled  with  non  collocated  transducers. 
Experimental  results  verified  the  non  self  adjointness  of  the  system 
due  to  non  collocated  control  and  a  procedure  was  developed  to 
identify  the  biorthogonal  eigenvectors  from  experimental  residue  data. 
Preliminary  verification  of  the  modal  post  processing  routine  was 
accomplished  using  experimental  data. 

A  model  of  a  collocated  piezoceramic  sensor  and  actuator 
occupying  a  subdomain  of  a  cantilever  was  developed  and  by  a 
minimization  of  a  selected  quadratic  functional  an  optimal  control  law 
was  derived.  This  control  law  was  evaluated  for  an  active  structure 
consisting  of  a  sensor  and  actuator  occupying  the  whole  domain  of  the 
beam  using  a  finite  difference  approach.  The  results  of  the  analysis 
indicated  that  a  step  Lyapunov  optimal  control  law  surpassed  a  rate 
control  law  in  reducing  the  response  of  the  beam  in  the  shortest  time. 
Also,  a  study  was  undertaken  to  evaluate  the  control  performance  of 
the  system  when  the  sensors  occupy  a  subdomain  of  the  structure. 
Results  of  this  study  indicated  that  through  the  sensor  output 
decreases  with  sensor  size,  control  can  be  regained  with  an  increase 
in  feedback  gain. 

Control  optimization  of  an  active  structure  using  a  limited 
output  state  feedback  control  was  investigated  using  piezoceramic 


sensors  and  actuators.  A  discrete  model  of  a  cantilever  structure 
with  multimode  control  was  formulated  in  state  space  and  the  output 
feedback  control  gain  matrix  optimized.  The  results  indicated  that  at 
the  optimum  gain  distribution  the  negative  real  part  of  the 
eigenvalues  is  increased. 

The  following  recommendations  are  presented  for  additional 
research  into  the  area  of  the  active  control  of  structures: 

a)  An  extension  of  the  finite  element  active  beam  element  to 
structures  of  higher  order  such  as  plates  and  shallow  shells  should  be 
undertaken.  This  will  permit  greater  flexibility  in  modelling 
structures  such  as  membrane  optics  or  large  dish  antenna  type 
structures. 

b)  Enhancement  of  the  dynamic  coupling  results  should  be  explored 
by  using  an  enhanced  structure  and  electronic  feedback  models  and 
identification  procedures. 

c)  Additional  experimental  work  should  be  performed  on  non 
collocated  control  from  both  a  system  stability  point  of  view  as  well 
as  identification  of  biorthogonal  modal  description .of  the  controlled 
response  of  simple  structures. 

d)  The  use  of  a  higher  order  control  compensator  incorporating 
narrow  band  filtering  characteristics  should  be  evaluated  for  optimum 
pole  placement  and  gain  distribution  for  multimode  control  of 
structure  using  a  limited  state  feedback  approach. 
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APPENDIX  A 


Piezoelectric  Constitutive  Equations 


Derivation  of  the  piezoceramic  constitutive  equations  has  been 
done  by  many  investigators®’®^’®^  ^  but  is  repeated  here  for 
completeness.  The  approach  commonly  used  and  followed  here  is: 

a)  Derivation  of  Poynting  equation  form  Maxwell's 

electromagnetic  theory  which  is  an  electromagnetic  statement  of  the 
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principle  of  the  conservation  of  energy 

b)  Using  a  modified  form  of  Poynting  equation  to  develop  the 
first  law  of  thermodynamics  for  a  piezoelectric  medium  ’ 

c)  Combining  this  law  with  crystallography  of  piezoceramic 

materials  to  construct  a  set  of  constitutive  equations  for 

8 

piezoceramic  transducers  ’ 


Derivation  of  Poynting  Equation. 


Maxwell's  equations  in  MKS  Units  are: 


7  X  H  =  D  +  3 

(Ampere's  Circuital  Law) 

(A-1) 

7  X  E  =  -B 

(Faraday's  Law) 

(A-2) 

O 

II 

too 

t> 

(Gauss's  Law  for  Magnetic  Field) 

(A-3) 

7D  =  pe 

(Gauss's  Law  for  Elect.  Field) 

(A-4) 

ig  =  charge  density 

[C/M] 

E  =  electric  field  intensity* 

[N/C], 

1 — 1 

> 

H  =  magnetic  field  intensity* 

[A/M] 

J  =  current  density 

[A/M^. 

[C/S  M^] 

D  =  electric  flux  density* 

[C/M^] 

B  =  magnetic  flux  density* 

[N  S/C  M] 

P  =  polarized  vector 

[C/M^] 

M  =  magnetized  vector 

[A/M] 

Eq  =  permittivity  of  a  vacuum 
Mg  =  magnetic  constant 

and  *  signifies  a  vector  field.  The  polarized  vector  is  a  measure  of 
the  electric  dipole  per  unit  volume  and  the  magnetized  vector  is  a 
measure  of  the  magnetic  dipole  per  unit  volume. 

Using  the  vector  identity 

7*  (AxB)=B-7xA-A*7xB  (A-7) 

and  equations  (A-1)  and  {A-2),  the  differential  form  of  Poynting's 
equation  is  derived  as  follows 


7  •  (E  X  H)  =  H  •  (7  X  E)  -  E  •  (7  X  H) 
7  •  (E  X  H)  =  -H  •  B  -  E  •  (j  X  5) 


(A-8) 

(A-9) 


Equation  (A-15)  is  the  integral  form  of  Poynting's  theorem  and  with 
units  [N  M/S]  =  [J/S]  =  [WATTS].  This  theorem  is  a  power  theorem 
representing  the  electromagnetic  statement  of  the  principle  of  the 
conservation  of  energy.  Each  part  of  equation  (A-15)  is  defined  as: 


Rate  of  change 
of  energy 
stored  in 
the  volume 


Net  inward  power  I  [Power  imparted  to 


flow  through  a 
closed  surface 


charges  in  the 
volume  enclosed  by  the 
surface  (Joule  heat) 


This  theorem  is  also  expressed  in  a  different  format  using  vector 
and  scalar  potentials  A  [Weber/Meter]  and  y  [Volt].  The  vector  field 
A  and  a  scalar  field  are  defined  from  the  knowledge 
of  J  and  then  the  field  vectors  0,E.B  and  H  are  computed  from  A  and  a. 
The  steps  necessary  to  obtain  Maxwell's  and  Poynting's  equations  in 
terms  of  electromagnetic  potentials  is  presented  in  reference  49. 

Expressed  in  terms  of  these  potentials  and  in  tensor 
notation 


{(£,6,  .H,B,)dV  =  jr,.  [r(D,  .  0,)  -  e,j^AH^]dS 
-fE,J,dV 


(A-16) 


with  the  Poynting  vector  defined  as 


Since  the  net  change  in  charge  within  the  body  of  an  insulator  is 
zero  (charge  resides  on  the  outside  of  the  body), 

Pe  =0 

Furthermore,  because  there  is  no  current  flow  in  an  insulator. 


Finally  since  the  materials  considered  in  this  thesis  are  polarized 
and  nonmagnetized  dielectrics 


M.  =  0 

From  the  above  conditions  Maxwell's  equations  for  the  piezoelectric 
material  are 


®ijk  ^k,j  ^  ^ 
®ijk  ^k,j  ®i 


Vs'.V 

\  s.* 

>  .*  V 


A 


(A-18) 

{A-19) 


Likewise  from  Poynting's  theorem  the  electromagnetic  energy 


flow  out  of  the  surface  enclosing  the  volume  of  interest  is 

Using  the  quasistatic  electric  approximation  Equation  (A-16) 
reduces  to 


-j  n.h.dS  =  |(E.D.  +  H.6.)dV 


{A-25) 


Using  the  principle  of  conservation  of  energy  and  equation  (A-25) 
the  equations  of  motion  for  a  piezoelectric  body  are  derived. 

For  a  piezoelectric  body  bounded  by  a  surface  S  with  a  unit 
outward  normal  n,  the  rate  of  increase  of  energy  (kinetic  and 
internal)  is  equal  to  the  work  done  by  the  surface  tractions  and 
the  flux  of  electric  energy  outward  across  S  assuming  no  body 
forces. 


3_  f  (1/2  u.u.  +U)dV  =  (f  (t.u.  -  n.yft.) 
Yt  Jy  J  J  J  J  J 


(A-26) 


where 

p  =  density  of  the  material 
u  =  internal  strain  energy 
tj  =  applied  tractions 

y  =  scalar  potential 

Dj  =  rate  of  charge  developed  per  unit  area 


Using  the  following  relations 


T..  .  =  ii.  Equilibrium 

I  J  9J  J 


D.  .  =  0  EquationA-21  for  the  dielectric  body 


E,.  =  -Y,i^  Definition 


S..  =  i{u,  .  +  u.  .)  Kinematics 

•w  IfJ  J*' 


where  T..  is  the  stress  tensor  and  S..  is  the  strain  tensor. 

I  J  I J 

O 

From  the  above  relations  and  equation  (A-26) 


0  =  T, .  5,J  -  £,6, 


(A-27) 


Equation  (A-28)  represents  the  first  law  of  thermodynamics  for  a 
piezoelectric  medium,  which  states  that  for  an  adiabatic  process,  the 
change  in  the  internal  energy  of  the  piezoelectric  material  is  equal 
to  the  rate  of  mechanical  and  electrical  work  done  per  unit  volume. 


Piezoceramic  Consititutive  Equations 


Starting  with  equation  (A-28)  an  electric  enthalpy  can  be 


defined  as 


;sa 


differentiating  equation  (A-28)  with  respect  to  time  and  combining  it 


with  equation  (A-27)  yields 


K  =  T,j  5,j  -  0,  E, 


(A-29) 


which  implies  that  H  is  a  function  of  the  strain  and  the  electric 
field.  Since 


a  _  afl  c 
H  -  -  S  -  * 

3S^.j  ij  3E.  'i 


subtracting  equations  {A-30)  from  (A-29)  results  in 


T  -  §  -  D  +  f  =  0 

^ij  °i  ^  aE.  ° 


S.j  and  E^  are  arbitrary  therefore 


T  -  lH_  n  -  ^ 

"ij  -  ’  °i  -  ■  aE. 


(A-30) 


{A-31) 


(A-32) 


for  equation  (A-32)  to  be  satisfied. 

Q 

Linear  piezoelectric  theory  ,  results  in  specifying  a  homogeneous 


quadratic  form  for  H; 


ijk  ^ij\  ■  ®ijK^i^jk  ’  ^  ^ij^i^j 


(A-33) 


where 


Cj.,.,  =  elastic  coefficients  for  a  constant  electric  field. 


=  piezoelectric  coefficients 


®i  jk 

e..  =  dielectric  coefficients 

'  \J 

From  equation  (A-32)  the  stress  and  the  electric  displacement  are 
specified  as 


T  C  .a  r  -  3H 

ij  "  ^ijkrkl  ®kij^k  “  3S.. 

°i  ®ikl  ^kl  ^ik  ^k  ^  IIt 


(A-34) 

(A-35) 


Substituting  equations  (A-34)  and  (A-28)  into  (A-33)  yields  in 


U  =  i 


^ijkl^ij^kl 


+  is 


(A-36) 


where  there  is  no  piezoelectric  interaction  term  in  the  positive 
definite  stored  energy  function  U. 

Equations  {A-33)  and  (A-34)  represent  the  constitutive  equations 
for  an  unbounded  piezoelectric  medium.  The  other  forms  noted  in  the 
literature  (6)  are: 
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f 

f 

) 

4 

I 


rn 


i 


I 


^  ^Ik°k 

^ij  ^  ^ijkl^kl  ■  ^kij°k 
"  '^ikl^kl  *  ^ik\ 


(A-40) 

{A-41) 

(A-42) 


Note,  the  superscripts  designate  the  following  condition: 

T  =  constant  stress 
E  =  constant  electric  field 
S  =  constant  strain 
0  =  constant  electric  displacement 

The  matrix  relationship  between  the  coefficients  in  equations 
(A-34)  to  (A-42)  are  expressed  in  the  following  matrix  from  as^®: 


[P^]  = 

[g]  =  CP^][cI] 

CS°]  =  [s^]  -  trans  ([d])[^'^][e] 
[e]  =  [d][c^] 

[h]  =  Cp^][e] 

[C°]  =  [C^]  +  trans  (Ce])[/J^][e] 
[P^]  =  [P^]  +  Ch][S°Jtrans  ([h]) 

=  [e^]  -  [d][C^]  trans  ([d]) 


I 

I 

i 

\ 

t 


To  use  equations  (A-33)  to  (A-42)  for  analysis  of  the 
electromechanical  coupling  of  unimorph  piezoceramic  sensors  and 
actuators  to  elastic  bodies  requires  knowledge  of  the  piezoceramic 
mechanical,  dielectric  and  piezoelectric  coefficients.  The 
piezoceramic  used  is  a  ferroelectric  material. 

To  obtain  the  necessary  coefficients,  for  the  particular 

g 

ferroelectric  used  requires  a  review  of  crystallography  and  its 
application  to  this  class  of  materials. 

Ferroelectric  ceramics  consist  primarily  of  a  single  crystal  type 
structure  of  polycrystalline  material.  When  this  material  is  poled  by 
the  application  of  an  intense  external  electric  field  at  elevated 
temperatures  it  is  made  piezoelectric  and  rendered  a  single  crystal 
symmetry  on  a  macroscopic  scale. 

Crystals  are  grouped  into  seven  systems  depending  on  their 
symmetry.  These  systems  are  further  divided  into  classes  according  to 
their  symmetry  to  a  point.  In  all,  there  are  32  classes  of  which  are 

Q 

considered  to  be  piezoelectric  . 

The  lead  zirconate  titanate  (PZT)  type  of  ferrelectric  is  a 
hexagonal  crystal  system  and  identified  using  the  international 
crystallographic  system  as  6mm  which  is  distinguished  by  an  axis  of 
sixfold  symmetry. 

To  relate  the  constitutive  parameters  to  the  different 

g 

crystalline  systems  an  elasto-electric  matrix  description  is  used  . 
This  method  of  description  relates  the  coupling  between  the  elastic 


and  dielectric  properties  (equation  A-33  and  A-34)  by  an  arrangement 
of  the  dielectric  and  elastic  coefficients.  The  elasto-electric 
matrix  consists  of  a  nine  by  nine  matrix  where  in  each  column  is 
related  to  the  independent  variable  such  as  stress  or  electric  field 
component  and  each  row  to  a  dependent  variable  such  as  strain  or 
electric  displacement. 

For  example,  the  triclinic  crystalline  system  or  the  one 
possessing  neither  symmetry  axes  or  symmetric  pianes  has  45  different 
coefficients  as  shown  below: 
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The  piezoelectricity  of  the  6inm  crystalline  material  is  now 
classed  by  three  independent  piezoelectric  constants,  e^^,  e^^  ^15 

or  as  noted  previously  d^^^,  d^^  and  d^g.  The  data  available  on  the 
piezoceramic  materials  is  usually  given  in  the  later  form. 

Equations  (A-33)  and  (A-34)  are  expanded  for  example  (ref.  74  )  to 
obtain  the  following  constitutive  equations  for  an  unbounded  6mm 
piezoceramic  medium. 


T, ,  =  C..U.  ,  +  C.nU/,...  +  C..,u.,  +  e.,,E.. 

11  11  1,1  12  2’2  13  3,3  31  3 


T«/,  =  C„^u,  ,  +  C«<u.,  +  C..,u.,  «  +  e.,,E., 
22  21  1,1  11  2,2  13  3,3  31  3 


^33  ~  Sl*^!,!  ^13^2,2  ^33“3,3  ®33^3 


^23  ■  ^44^‘^3,2  *  *^2,3^  ®15^2 


^31  ‘  ^44^‘^3,1  ‘^1,3^  ®15^1 


^12  ■  ^66^*^!, 2  ‘^2,1^ 


°1  ^  ®15“3,1  ®15‘'l,3  ■  ^11^1 


°2  ■  ®15^“3,2  “2,3^  '  '^11^2 


°3  ■  ®3l‘^l,l  ®3l“2,2  ®23“3,3  '  ^33^3 


List  of  Piezoelectric  Coefficients  and 

Their  Units 

ijkl.^pq 

Elastic  Stiffness  Constant 

N/M^ 

ijk/ip 

Piezoelectric  Constant 

M/V  or  C/N 

ijk,®ip 

Piezoelectric  Constant 

C/M^ 

i 

Electric  Field  Component 

V/M 

ijk.^ip 

Piezoelectric  Constant 

V-M/N  or  M^/C 

ijk.^ip 

Elect.  Enthalpy  Density 

J/M^ 

i 

Electric  Disp.  Component 

C/M^ 

Capacitance 

Farad 

ijkl.^pq 

Elastic  Compliance  Constant 

M^/N 

ij 

Impermittivity  Constant 

M/Farad 

0 

Perm,  of  free  space  (8.85e“12) 

Farad/M 

1  1 

Perm,  component 

Farad/M^ 

aaaaaa  aacgsass?  >j^xcvi  ^ 
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APPENDIX  B 


FINITE  ELEMENT  ACTIVE  BEAM  MODEL 


This  appendix  contains  the  following: 

a)  The  derivation  of  the  element  matrices  of  a  beam  finite 
element  model  with  a  piezoceramic  driver. 

b)  The  derivation  of  the  element  matrices  of  a  beam  element 
with  a  piezoceramic  sensor  and  actuator. 

c)  Investigation  of  the  effects  of  Guyan  condensation  on  the 
active  structural  damping  matrix. 


BEAM  ELEMENT  WITH  AN  ACTUATOR 


The  finite  element  model  °  of  a  beam  element  shown  in  figure  B-1 
with  a  perfectly  bonded  unimorph  piezoceramic  driver  and  subjected 
only  to  bending  is  derived  follows:  Fixing  node  i-1. 


“1-1  =  ®i-i  = 


Equilibrium  for  the  beam  is 


^  =  It  C-M,  *  V,(L-x)  *  K(,Vp] 


where  Vq  is  the  DC  voltage  applied  to  the  driver. 


(B-1) 


•  ■  •  «  •  «  •  •  •  •  •  »  '  m  •  ■  *  *  •*  »**^^**  '  •*  •  •  • 


Actuator  Element  With  Free  Body  Diagrams 
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Integrating  once  gives 


H  =  Ei  r>  *  "d  * 'l 


(B-2) 


x=0,  ^=0 


X  =  L  ^ 
^  dx 


■®i 


Solving  for  the  rotation  at  node  i  gives 


V.L" 


e,  =  ^  [-n,L  -  4-  -  Vo'-] 


(B-3) 


Integrating  equation  (B-2)  gives 


x=0,w=0  C2=0 


2  ‘'2 


] 


(B-4) 


X  =  L,  w  =  w. 


Solving  for  the  displacement  at  node  i  gives 


v!.  w!l 

2  3  2  J 


(B-5) 


Combining  equations  (B-3)  and  (B-5)  in  matrix  form  and  taking  the 


inverse  of  the  flexibility  matrix  gives 


Fixing  node  i  results  in  the  following  free  body  diagram  shown  in 
figure  B-2.  Equilibrium  for  the  beam  is  now  with  W.=  0  and  0.=  0. 


[Vi.ix 


^-1 


^  W 


(B- 


where  Vq  is  again  the  DC  voltage  applied  to  the  driver.  Integrating 
once  gives 


dx  ■  El  [”(-1*  *  'z'  *  Vo*  * 

x  =  L.f  xO  /.c, 


Therefore  solving  for  the  rotation  at  node  i-1  gives 


located  Rate  Control  Example  Configuration 


<<s<: 


-  Kpl/pj/L  )  L  -1  -1  J  <  (M._^  +  KjjVp)/L  j 


Substituting  equation  (B-12)  into  (B-15)  gives 


-  W/'- 


«i-1 


0i.ll- 


{B-16) 


where 


Combining  equations  (B-6),  B-12),  (B-14)  and  (B-16)  gives  the  final 
result 


''i-1 

'  ^-1 

12 

{«i-i  ^  w^'- 

>  _  El 

-  6 

4 

Sym 

' 

®i-lL 

''i 

-12 

6 

12 

”i 

-  6 

2 

6 

4  _ 

^  (Ml  -  K,Vo)/L  ^ 

®i‘- 

N  > 

(B-17) 


^  V 


as  in  section  2.2.  The  free  body  diagrams  in  figure  B1  and  B2  are 
val id. 


From  equation  (2.22)  in  chapter  II,  the  equation  for  the  voltage 
form  the  sensor  is 


V.{t)  =  K  f  ^  dx 

at 


where 


(B-19) 


From  chapter  11,  rate  feedback  is  defined  as 


(I 


Combining  equation  (B-20)  with  the  integral  (B-19)  evaluated  and 
substituting  for  in  equation  (B-18)  yields 


Expanding  the  right  hand  side  of  equation  (B-21)  gives  the  final  result 


■=  H[K]  +  KK 

0  0  0  0 

0  10-1 

0  0  0  0 

- 

r  •  1 

^-1 

^-1 

L 

0-101 

(B-22) 


GUYAN  CONDENSATION  OF  THE  ACTIVE  DAMPIWS  MATRIX 


Identification  of  the  dynamic  coupling  of  the  active  was 
presented  in  section  chapter  III.  The  five  degree  of  freedom  model 
consisted  of  only  translational  degrees  of  freedom  with  the  rotational 
d.o.f.  condensed  out  using  Guyan  condensation.  Presented  in  this 
section  is  the  results  of  condensation  of  the  small  four  element  of  a 
beam  shown  on  figure  B-3  with  two  collocated  sensor  actuator  pairs. 

The  structure  considered  was  assumed  not  to  have  any  internal 
damping.  The  mass  ans  stiffness  properties  of  the  structure  can  be 
arranged  as 


Given  equation  {B-23)  a  transformation  matrix®®’®^  can  be  developed 
such  that 
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(B-24) 


where 

[T]  =  -  (B-25) 

Using  equation  (B-25),  equation  (B-23)  can  be  transformed  to 
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where 

and  [K]=  [[I]  [T-f  IwHlj 


CM]  =  [I]  [T]^  CM] 


Cl] 


With  equation  (B-28)  partitioned  as  shown  in  equation  (B-23)  the 
transformation  matrix  can  be  computed  according  to  equation  (B-25) 
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1.6080 

-1.2680 

To  condense  the  structural  electronic  damping  matrix  to  only 
displacement  d.o.f.  requires  the  following  operation. 


[C]  =  [ci]  ct]^][Cep]  j'cnj 


The  result  of  this  operation  is 


(B-29) 
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(B-31) 

The  condensed  form  of  the  electronic  damping  matrix  is  fully  populated 
and  symmetric  for  collocated  sensors  and  drivers. 


APPENDIX  C: 


System  Identification  of  Vibrating  Structures 

System  identification  is  the  process  of  using  measured  test  data 
of  a  structure's  response  and  usually  some  a  priori  information  of  the 
system's  parameters  to  develop  a  mathematical  model  which  simulates 
the  system's  behavior.  System  identification  can  be  accomplished  in 
the  time  or  frequency  domain  employing  either  physical  or  modal  data. 
Common  approaches  used  are: 

a)  The  minimization  of  the  error  satisfying  the  system's 
equations  with  the  measured  response. 

b)  The  minimization  of  the  differences  between  the  measured  and 
analytical  response. 

c)  Seeking  a  new  set  of  system  parameters  that  differ  the  least 
from  a  priori  or  specified  values  using  the  system  equations  as 
constraints. 

Given  the  matrix  description  of  a  system  as 

[M]  {X.}  +  [C]  {X.}  +  [K]  {X.}  =  {f.}  (C-1)  ‘ 
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where  the  matrix  coefficients  have  been  identified  and  represents  the 
vector  of  the  physical  degrees  of  freedom  and  {f^.}is  the  vector  of  the 
forcing  functions,  equation! C-1 ) becomes  the  baseline  analytical  model. 
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This  model  can  be  used  for  example  to  study  changes  in  the  physical 
system  as  well  as  changes  in  the  boundary  conditions. 

The  system  can  also  be  characterized  by 

[Z(a.)]  {x.{u,)}  =  {f.(u.)}  (C-2) 

where  the  inverse  of 

[Z(u.)]‘^  =  [-a)2[M]  +  MC]  +  [K]‘^]=  [Y(a.)]  (C-3) 

is  the  mobility  matrix  and  each  element  is  a  transfer  function.  The 

identification  of[Y(a>)]  can  be  used  to  study  the  effects  of  various 

65 

forces  on  the  structure  or  for  a  component  synthesis  analysis  . 

The  system  can  also  be  characterized  by 

A?  [M]  {9^}  +  A^CC]  {♦.}  +  [K]  =  {0}  (C-4) 

where  and  are  the  eigenpararoeters.  In  this  form  the  general 
validation  of  an  existing  analytical  model  can  be  obtained.  This 
particular  characterization  was  used  in  an  algorithm  developed  to 
identify  the  baseline  system  matrices  without  active  control.  A  brief 
outline  of  this  procedure  is  given  next  for  completeness. 

The  system  matrices  of  a  dynamic  system,  such  as  the  one 
described  by  equation  (C-1),  satisfy  the  eigenvalue  problem 

A?  K*,  =  0  (  i  =  1,2 _ 2n)  (C-5) 


where 


X.  =  1th  complex  eigenvalue 
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If  one  uses  arbitrary  matrices  in  place  of  the  correct  system 
matrices  in  Equation  (C-5),  it  is  not  satisfied  exactly.  The  norm  of 
the  resulting  error  matrix  is  then  minimized  over  all  the  modes  with 
respect  to  the  matrix  elements  subject  to  symmetry  constraints.  This 
process  results  in  a  set  of  homogeneous  simultaneous  equations  given 


MLj  +  LjM  +  CL2  +  L2^C  +  KL3  +  L3^K  =  0 
ML2T  +  L2M  +  CL^  +  L^C  +  KLg  +  =  0 

ML3T  +  L3M  +  CLg^  +  LgC  +  KLg  +  Lg  K  =  0 


where 


L  =  ZF^(F^)^  Z^, 
L  =  Z(F^)2z, 

L  =  ZfV, 


z  =  CRe(*.) 
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Ar  =  cliag[Re(«. )  ...  Re(4.)] 

Aj  =  diag[Itn(*. )  ...  Im(*. )] 

If  some  of  the  matrix  coefficients  are  know  a  priori,  the 
remaining  unknown  coefficients  can  be  determined  from  these  equations. 

If  the  identified  matrices  are  such  that  they  satisfy  the 
eigenvalue  problem  with  respect  to  the  measured  modal  parameters 
exactly,  the  orthogonality  conditions  are  also  satisfied.  If  the 
error  in  the  eigenvalue  problem  is  not  zero,  orthogonality  is 
satisfied  approximately.  However,  even  in  this  case,  orthogonality  is 
satisfied  exactly  with  respect  to  the  identified  modal  parameters, 
which  in  turn  are  approximately  equal  to  the  measured  modal 
parameters. 

The  identification  procedure  used  to  compute  the  mass,  damping 
and  stiffness  matrices  required  that  some  of  the  mass  or  stiffness 
elements  to  be  known  a  priori.  To  provide  this  information,  models  fo 
the  various  test  cantilevered  beams  were  developed  using  the 
GTSTRUOL®^  finite  element  code. 

The  finite  element  model  consisted  of  plane  frame  XZ  type 
elements.  The  mass  of  the  accelerometer  located  at  the  tip  of  the 
test  beams  was  included  as  a  lumped  mass  at  the  tip  node  point  of  the 
model.  The  mass  of  the  piezoceramic  transducers  and  wires  was 
considered  negligible  and  therefore  not  included.  Finally,  depending 


APPENDIX  D: 


ACTIVE  CONTROL  EXPERIMENTS 


Experimental  Setup  and  Equipment 

Structural  dynamic  control  experiments  were  undertaken  to 
evaluate  the  performance  of  the  active  control  of  simple  structures. 
Three  cantilever  beam  were  selected  which  had  the  following 
dimensions: 

#1)  Steel  Beam  24x2x.5  inch. 

#2)  Aluminum  Beam  24x2x.25  inch. 

#3)  Aluminum  Beam  9x.651x.l72  inch 

Beam  #3  was  used  for  the  identification  of  the  dynamic  coupling 
coefficients.  All  three  beams  were  welded  to  a  thick  metal  plate 
which  was  bolted  a  support  frame.  Figure  D-1  is  a  drawing  of  a 
typical  beam  experimental  setup. 

Piezoceramin  unimorph  transducers  were  used  to  sense  and 
dynamically  control  the  cantilevered  beams.  The  transducer  locations 
for  the  three  beams  are  shown  in  figure  D-2.  The  transducers  were 
supplied  by  Piezoelectric  Corp.  and  were  made  out  of  G1195  material 
which  is  a  lead  zirconate  titanate  (PZT)  type  compound.  Depending  on 
the  transducer,  it  consisted  of  either  3/4  x  3/4  x  .010  in.  or  a  1.5  x 
.655  X  .010  thin  sheet  of  piezoceramic  sandwiched  by  two  nickel 


9.625 
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electrodes  that  cover  the  entire  surfaces.  Shielded  leads  were 
soldered  to  the  electrodes  at  the  corners  of  the  transducers.  Table 
D-1  lists  the  material  and  piezoelectric  specification  for  the  PZT 
material . 

For  each  configuration,  the  transducers  were  bonded  to  the  beam  a  at  the 
desired  locations  with  a  thin  layer  of  epoxy..  A  small  amount  of  micro 
bubbles  were  mixed  in  the  epoxy  to  keep  the  bottom  of  the  electrodes  were 

insulated  from  the  beam  surface  and  to  maintain  an  even  layer  of  adhesive. 

min. 

The  epoxy  consisted  of  a  two  part  mix  which  required  fiveHo  cure  at  room 
temperature. 

In  general,  the  response  measurements  of  the  beams  were  made 
using  accelerometers  (BBN  model  507  Sn  1346,  PCB  model  303A  SN 
7384).  The  accelerometers  were  usually  fixed  at  a  single  location 
(tip)  and  an  impact  hammer  with  a  force  transducer  (PCB  model  208A03 
SN  4417  and  PCB  086B80  Snl284)  was  moved  along  the  beam  to  the 
different  locations  to  apply  impact  excitation. 

For  each  collocated  transducer  pair,  one  transducer  was  selected 
as  the  sensor  and  the  other  as  an  actuator.  The  current  output  of  the 
sensor  was  conditioned  by  an  active  feedback  control  system  which 
provided  high  voltage  output  to  the  actuator. 

Two  prototype  active  feedback  control  units  were  provided  by  Dr. 

R.L.  Forward  from  Hughes  Research,  Malibu  CA.  These  units  were 
designed  for  the  sole  purpose  of  demonstrating  the  potential 
performance  of  active  control  using  piezoceramic  transducers.  The 
units  basically  consist  of  a  preamplifier,  filter  and  high  power 


Table  D-1.  G-1195  Piezoceramic  Coeff icienba . 


Electrodes  Nickel 
Modulus  =  6.3E*10  N/M^ 

Capacitance  (.75  x  .75  x  .01  in.>  =  2.2E-08  C/V 

Capacitance  <1.555  x  .6036  x  .01  in.)  =  3.9E-08  C/V 

931  =  1.14E-02  VM/N 

d3i  =  1.79E-10  M/V 

€33  =  1.65E-0a  Farad/M 

Curie  Temp.  =  360*  C 

Density  =  7.6  gr-Zcc 
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output  sections.  The  preamplifier  section  is  simply  a  charge 

amplifier  which  converts  the  output  of  a  piezoceramic  thin  sheet 

ceramic  strain  gage  (current)  to  a  voltage.  The  voltage  is  then 

filtered  to  a  desired  frequency,  phase  and  frequency  bandwidth  (Q) 

prior  to  being  amplified  in  the  high  power  section  for  driving  other 

piezoceramic  strain  gages.  Figure  D-3  is  a  drawing  of  the  front  panel 

of  an  active  control  unit.  A  detailed  analysis  of  the  active  feedback  and 

their  calibration  is  covered  later. 

To  measure  and  reduce  the  dynamic  response  data  a  16  channel 

Computer  Aided  Test  System  (Genrad  Model  2515  was  used).  This  system 

provided  a  multitude  of  methods  to  evaluate  the  quality  of  the  data  as 

well  as  computing  the  change  in  modal  damping  upon  the  application  of 
51 

active  control  .  Figure  0-4  shows  a  typical  experimental  setup  with 
all  of  the  primary  signal  conditioning  equipment. 


Test  Procedure 


Prior  to  each  test,  the  gain,  center  frequency,  bandwidth  and 
phase  of  the  active  damping  units  was  adjusted.  The  Q  or  bandwidth  of 
the  filter  circuit  was  set  at  its  narrowest  setting  which  was  20. 

Using  the  output  of  the  accelerometer  or  a  sensor,  the  frequency  at  a 

particular  gain  setting  was  fine  tuned  until  a  minimum  output  was 
obtained  for  the  given  sinusoidal  excitation  applied  by  an 
piezoceramic  actuator.  The  optimum  phase  setting  was  at  180  degrees 

for  collocated  sensors  and  was  not  optimized  for  the  non  collocated 

test. 
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After  the  control  units  were  set  for  a  particular  test 

configuration,  impact  response  tests  were  undertaken.  Both  baseband 
51 

and  zoom  frequency  response  analysis  were  used  depending  on  the 
objective  of  the  particular  test. 

The  impact  and  response  signals  obtained  for  zoom  analysis 
required  filtering  which  was  done  by  a  pair  of  Krohn-Hite  dual  channel 
filters  (models  3343  and  3342R).  The  filtering  was  necessary  to 
maintain  an  acceptable  signal  to  noise  ratio.  A  wideband  random 
signal  was  simultaneous  applied  to  both  channels  and  a  transfer 
function  was  computed  between  them  to  measure  any  phase  shift  caused 
by  the  filters.  No  phase  shift  was  detected  in  the  frequency  band  of 
the  experiments. 


DATA  PROCESSING  AND  REDUCTION 


Using  transfer  function  analysis  techniques  the  modal  parameters 

for  the  beams  were  determined.  The  impact/response  measurements  were 

averaged  together  together  for  each  beam  location  and  stored  for  later 

analysis.  When  a  particular  test  condition  was  completed,  the 

resulting  transfer  functions  were  processed  using  SDRC  Modal  Plus 
51 

software  to  obtain  complex  modal  parameters  and  mode  shapes. 

Various  options  available  for  modal  analysis  were  explored  and 
compared.  Differences  between  the  parameters  identified  on  this 
single  input  output  type  system  were  slight  and  it  was  finally  decided 
that  the  complex  exponential  technique  was  adequate  to  obtain  the 
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complex  pole  and  residue  data.  Finally,  a  circle  fit  routine  was  used 
to  obtain  the  real  and  imaginary  modal  coefficients. 

MODEL  OF  THE  FEEDBACK  ACTIVE  CONTROL  ELECTRONICS 


Schematic  drawings  were  supplied  with  the  prototype  electronics 
and  are  shown  in  figure  D-5.  Two  tasks  were  undertaken  to  model  the 
control  units.  The  first  task  was  to  perform  a  detailed  system 
analysis  using  the  schematics.  Secondly,  since  the  possibility 
existed  that  components  were  altered  from  the  original  design,  a 
series  of  calibratioin  tests  as  well  as  physical  examination  of  the 
internal  components  were  performed.  The  test  also  served  to  further 
verify  the  differential  equation  model  of  the  active  control  feedback 
filters. 

As  mentioned  previously,  each  active  control  feedback  filter 
consisted  of  a  zero  impedance  preamplifier,  bandpass  filter,  phase 
shifting  circuit  and  high  voltage  output  amplifier.  Each  subsection 
is  identified  on  the  schematic  drawing  in  figure  D-5. 

The  preamplifier  subcircuit  shown  in  figure  D-6  can  be 
apporpriately  modelled  as  shown  in  figure  D-7.  The  input  resistance 
is  selectable  from  10®  to  10^^  ohm  with  gain  knob  #1.  The  corner 
frequencies  of  the  various  parallel  capacitor  resistor  networks  were 
not  within  the  30  to  600  HZ  experimental  operating  range  and  therefore 
their  effect  was  considered  negligible. 


The  equation  modelling  this  portion  of  the  control  unit  is 


%(<-) 

The  bandpass  filter  section  is  drawn  in  figure  D-8.  The  filter 
center  frequency  and  Q  were  adjustable.  The  filter  Q  was  set  at  its 
maximum  value  of  20  for  all  experimental  tests.  Figure  D-9  is  a 
impedance  block  drawing  of  figure  D-8  to  assist  in  evaluating  this 
part  of  the  control  unit.  Using  the  principle  of  superposition,  the 
solution  for  the  voltage  at  node  3  is 

V.(a.)  =  -  V.  V,  -  V.  (D-2) 

Z,  Z5(Xj)  Xg 

where  x^  represents  the  position  of  the  Q  adjustment  knob.  Solving 
for  the  voltage  at  nodes  5,  7  and  8  yields 


V5  =  -(Z^/Z3(X2))V3 

(D-3) 

V7  =  -(Z/Zg(x2))V5 

(D-4) 

1/3  = 

(D-5) 

where  X2  represents  the  position  of  the  center  frequency  adjustment 
knob. 

Since  Zq  =  Z^in  equation  (0-5) 
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Combining  equations  (D-3),  (D-4)  and  (D-5)  with  (D-2)  yields 


Vg  (o))  =  0.1a 
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b  “  0)  +  jua 
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=  V6p{a.) 
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a  =  oj  /Q,  b  =  0) 
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Q  =  Zg/22(Xj) 

The  output  from  the  bandpass  filter  is  conditioned  by  the  phase 
shifting  network  shown  in  figure  D-10.  The  phase  did  not  need 
adjustment  a  collocated  sensor  and  actuator.  The  phase  was  set  at 
zero  and  toggle  switch  which  allowed  a  180  degree  shift  was  used  to 
set  the  proper  polarity  for  active  feedback  control. 

The  transfer  function  of  the  phase  shifting  network  between  nodes 
5  and  10  is 

VjQ(a»)  =  -  [1/(1  +  j(9.4  X  lO'^a,)]  (D-8) 


The  gain  was  computed  for  the  transfer  function  in  equation  (D-8)  for 
a  frequency  range  of  1  to  5000  rad/sec.  The  results  are  plotted  in 
figure  D-11  and  indicate  a  rolloff  at  about  2000  rad/sec. 


After  the  phase  network,  the  signal  is  further  conditioned  by  a 
high  gain  voltage  amplifier  network  shown  in  figure  D-12.  Analysis  of 
this  portion  yielded 
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The  parameter  a^  varies  between  1.352x10^^  to  1.375x10^  as  the  gain 
adjustment  know  #2  is  adjusted  from  1  to  10.  Using  the  extreme  values 
of  a^  the  gain  was  computed  for  a  frequency  range  of  1  to  5000 
rad/sec. 

These  results  were  presented  in  figures  D-13  ,  and  indicate  a 

gain  range  of  .3  to  3.3  through  this  section  the  final  part  of  the 
amplifier  had  a  gain  of  30. 

Combining  equations  (D-8),  (D-9)  and  (D-10)  and  ignoring  the 
poles  and  zeros  out  of  the  range  of  operation  yields 


(D-10) 


where 

G  is  the  total  gain  through  the  system 


Vj(u))  =  "Rpii|y|(w) 


The  range  of  the  voltage  gain  from  equation  (D-8),  (D-9)  and 
(D-10)  is  approximately  .9  to  9.9  when  operating  between  10  to  2000 
rad/sec.  To  actually  determine  the  gain  for  a  control  experiment 
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required  that  each  control  unit  be  calibrated  after  a  particular  test. 
Calibration  was  accomplished  by  providing  a  wideband  current  input  to 
the  preamplifier  using  a  resistor  in  series  with  random  voltage 
source.  The  transfer  function  between  the  preamplifier  output  voltage 
and  the  high  voltage  output  provided  the  data  for  a  least  squares 
identification  routine  which  identified  the  gain,  Q  and  center 
frequency.  This  technique  was  presented  in  chapter  III. 

Table  D-3  is  the  peak  voltage  at  the  center  tuned  frequency  when 
a  1  random  signal  is  passed  through  a  1  x  10®  ohm  resistor  to  the 


preamp  input. 


EXPERIMENTAL  RESULTS 


A  variety  of  different  tests  were  done  to  examine  the  active 


electronic  damping  performance  using  the  G-1195  PZT  type 
sensors/drivers  on  the  cantilever  beams  and  Table  D-2  lists  the  major 
tests  conducted. 

Tests  1-23  were  exploratory  in  nature  and  will  not  be  discussed 
further.  Test  25-26  were  active  control  experiments  applied  to  beam 
#3  and  were  used  for  the  identification  of  the  dynamic  coupling 
coefficients. 

Prior  to  each  active  control  test,  the  baseline  (no  control) 
eigenvalues  and  eigenvectors  were  obtained  for  each  configuration. 


Figures  D-14a  &  D-14b  are  a  typical  sample  of  the  transfer  function 
data  of  the  first  and  second  modes.  Table  D-4  lists  the  complex 
eigenvalues  for  the  first  ten  modes  of  beam  #3  and  the  experimental 
eigenvectors  for  the  first  two  bending  modes.  The  eigenvectors  were 
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Table  D-2.  Active  Control  Experiments 


Test  Mode 


Date  Remarks 


Bean  #1 

llDecd4 

No  Active 

Control 

llDecd4 

S#2y'A#2 

llDec84 

S#2/A#2 

llDecS4 

S#i/A#l 

12Dec&4 

S#l/A#l 

13DecB4 

S#l/A#l  & 

S#2/A#2 

&Mar85 

No  Active 

Control 

8Mar85 

S#2/A#2 

6Mar85 

S#2/A#2 

7Mar85 

S#2/A#2 

8Mar85 

S#2/A#2 

8Mar85 

S«2/A#2 

Jun85 

No  Active 

Control 

Jun85 
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Figure  D-15b  Circle  Fit  Second  Mode  -  No  Active  Control 


Figure  D-16  Baseline,  Low  And  High  Gain  Comparisons 
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